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CHAPTER  I 


INTRODUCTION 

Electron  acoustic  grades  have  been  used  in  experimentally  deter¬ 
mining  the  electron  density  in  the  plasma  sheath  surrounding  a  reentry 
vehicle  (Ref.  1).  However,  these  methods  employed  resonance  effects 
which  are  in  observable  in  plasmas  with  large  electron-neutral  particle 
collision  frequencies.  A  short-pulse  diagnostic  technique  using 
electroacoustic  resonances  and  time  domain  analysis  of  the  reflected 
pulse  i«  described  by  Lustlg,  Baird  and  Evald  (Ref.  2).  For  slender 
conical  reentry  vehicles  which  reenter  the  earth's  atmosphere  at 
hypersonic  velocities,  the  plasmas  of  interest  are  collision-dominated. 
When  the  vehicle  readies  lower  altitudes,  the  collision  frequency 
can  be  larger  than  the  operating  frequency  of  the  electron  acoustic 
probe.  The  plasma  electron  density  is  inhomogeneous  and  Increases 
from  some  value  at  the  skin  of  the  vehicle  to  a  peak  value  at  a 
distance  typically  1  an  from  the  skin.  Rather  than  using  resonance 
effects,  it  is  proposed  to  use  propagating  electron  acoustic  waves 
to  probe  this  region.  The  waves  can  be  generated  at  the  boundary 
between  the  vehicle  and  the  plasma  by  an  electromagnetic  wave  incident 
obliquely  upon  the  boindary  from  a  dielectric  material.  In  the 
inhomogeneous  plasma  region  the  acoustic  waves  are  reflected,  due  to 
the  electron  density  gradient,  and  return  to  the  boisidary.  The 
reflected  acoustic  wave  which  returns  to  the  boundary  will  be  con¬ 
verted  into  a  transverse  electromagnetic  wave  which  can  be  detected 
(if  it  has  sufficient  saqilitude).  The  acoustic  wave  travels  much 
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slower  than  an  electromagnetic  wave  and  the  returning  pulse  of  acoustic 
energy  will  lag  the  returning  electromagnetic  pulse  In  time.  It  is 
proposed  to  determine  the  electron  density  profile  from  Its  effect  on 
the  acoustic  waves  by  observing  the  reflected  acoustic  waves. 

It  is  shown  later  (Fig.  22)  that  in  a  homogeneous  plasma  which 
is  collision  dominated  the  acoustic  waves  are  attenuated  very  rapidly 
and  are  not  likely  to  be  observable  for  propagation  distances  larger 
than  0.05  mm  at  the  most.  In  an  inhomogeneous  plasma,  the  acoustic 
waves  are  coupled  to  the  electromagnetic  wave  and  the  acoustic  waves 
may  not  be  so  severely  attenuated.  This  investigation  is  concerned 
with  a  steady-state,  plane  wave  analysis  of  tibia  coupling  effect  and 
its  implications  on  electron  acoustic  wave  propagation  in  an  inhomoge¬ 
neous  plasma. 

Electron  acoustic  waves  depend  upon  the  dynamic  properties  of  the 
plasma.  For  plasmas  with  sufficiently  large  electron  densities,  it  is 
appropriate  (Ref.  3)  to  describe  the  acoustic  waves  in  the  plasma  by 
a  scalar  pressure  field.  The  scalar  pressure  field  can  describe  a 
propagating  pressure  wave  and  is  physically  determined  by  the  dynamic 
properties  of  the  plasma  and  the  source  exciting  the  wave. 

Acoustic  wave  propagation  in  a  plasma  is  possible  only  at  a  finite 
plasma  temperature.  A  cold  plasma  (Incompressible  plasma)  is  one  with 
zero  temperature,  and  no  motion  of  the  particles  in  the  plasma  is 
assumed.  A  warm  plasma  (compressible  plasma)  exhibits  a  finite 
temperature  with  corresponding  random  thermal  motijns  of  the  particles 
in  the  plasma.  A  hot  p.’asma  is  one  with  extremely  high  temperatures 


(e.g.,  full on  temperatures) .  Electron  acouetic  wave  propagation  will 
be  considered  for  the  wan  plaaaa  only.  The  propagation  of  an  acoustic 
wave  froa  oi.e  region  of  the  plaaaa  to  another  region  ia  a  consequence 
of  the  raodoa  thermal  notions  of  the  electrons.  This  makes  the  local¬ 
isation  of  a  disturbance  impossible,  and  the  electrons  carry  the  dis¬ 
turbance  froa  one  region  to  another  (Sef.  3).  The  average  effect  is 
not  Isotropic  even  though  the  notions  of  the  electrons  are  random. 

The  electrons  moving  in  the  direction  of  the  wave  will  experience  a 
larger  change  of  momentum  and  there  will  be  a  net  tendency  to  carr- 
the  disturbance  in  the  direction  of  the  wave  (Ref.  3).  It  is  important 
to  distinguish  between  the  random  motions  of  the  individual  electrons 
and  the  collective  oscillatory  motion  of  the  electron  plasma  considered 
as  a  me  dim.  The  me  dim,  in  this  sense,  ia  described  as  a  single- 
fluid  electron  gas.  The  dynaalc  effects  give  rise  to  the  collective 
motion  of  the  plasma  and  form  the  physical  basis  for  acoustic  wave 
propagation.  The  collective  motion  of  the  plasma  medium  is  described 
in  a  macroscopic  manner  by  a  hydrodynamic  formulation  of  the  equations 
describing  wave  propagation  in  a  plasma. 

TWo  approaches  used  in  solving  plasma  problems  are:  a  micro¬ 
scopic  gas-kinetic  treatment  using  the  Boltzmann  transport  equation 
together  with  Maxwell's  equations  of  electrodynamics;  and  a  macro¬ 
scopic,  hydrodynamic  treatment  using  the  equations  of  conservation  of 
mass  and  momentum  together  with  Maxwell's  equations  (Ref.  4).  The 
kinetic  treatment  is  very  difficult  to  use  in  mathematically  modelling 
the  modes  of  plasma  oscillations  without  the  imposition  of  serious 
physical  restrictions  In  order  to  simplify  the  mathematics  in  the 


kinetic  treatment,  it  is  necessary  to  make  asswptions  which  would  make 
it  more  reasonable  to  use  the  simpler  hydrodynamic  treatment  (Ref.  A). 
The  fimdamental  equations  of  hydrodynamics  can  be  obtained  from  the 
Boltzmann  transport  equation.  To  a  certain  extent,  then,  the  hydro- 
dynastic  treatment  is  an  approach  to  the  correct  kinetic  treatment 
(Ref.  4).  It  must  be  stated,  therefore,  that  the  hydrodynamic  treat¬ 
ment  used  in  this  study  is  only  an  approximation  to  an  exact  solution 
of  the  present  plasma  problem.  It  mhould  be  noted,  however,  that  the 
hydrodynamic  approach  fails  intrinsically  only  for  cases  requiring 
more  than  a  phenomenological  approach  to  damping  effects  caused  by 
collisions  of  the  plasma  particles  and  for  problems  in  which  the 
velocity  distribution  function  is  specifically  involved  (Ref.  4). 

Such  cases  might  Include  calculations  of  scattering  cross  sections  or 
of  the  thermal  conductivity  of  a  plasma.  In  the  present  problem,  we 
are  concerned  with  macroscopic,  organized  behavior  of  a  plasma  median 
in  which  the  collision  frequency  is  energy-independent  and  is  described 
in  a  phenomenological  manner.  The  effect  of  the  collision  frequency 
is  included  as  a  drag  force  on  the  ordered  motion  of  the  electrons  in 
the  conservation  of  m omen  tun  equation.  The  hydrodynamic  treatment  in 
this  case  should  be  an  entirely  justified  approach.  There  are,  however 
certain  effects,  such  as  Landau  damping,  w  lich  are  not  included  in 
the  hydrodynamic  treatment.  These  will  be  discussed  in  detail  in 
Chapter  II. 

Electron  acoustic  waves  can  be  excited  by  an  electromagnetic 
wave  incident  obliquely  upon  a  plane  boundary  between  a  uniform 
dielectric  material  and  a  plasma.  Since  the  boundary  conditions  are 


both  electromagnetic  and  acoustic,  the  incident  electromagnetic  wave 
will  excite  an  acoustic  wave  at  the  botadary.  It  is  necessary,  however, 
that  tiie  electromagnetic  wave  have  a  component  of  the  electric  field 
strength  vector  perpendicular  to  the  boundary  (Ref.  5).  For  the  two- 
dieensional  geometry  considered,  this  means  that  the  incident  electro¬ 
magnetic  wave  nay  be  taken  to  be  vertically  polarized  with  the  electric 
field  strength  vector  in  die  plane  of  Incidence  in  the  most  general 
case  of  interest.  This  is  discussed  further  in  Chapter  III. 

A  survey  of  the  literature  was  made  to  determine  as  completely  as 
possible  the  previous  work  done  in  die  area  of  wave  propagation  in 
compressible  plasmas.  A  significant  result  was  the  fact  that  very 
few  mmterlcal  results  have  been  ptAlished  which  describe  the  propaga¬ 
tion  of  an  electron  acoustic  wave  in  an  Inhomogeneous,  lossy  and 
botmded  plasma.  The  few  numerical  results  published  were  obtained 
for  wave  propagation  in  the  ionosphere  (Ref.  6).  The  plasma  under 
consideration  in  this  report  is  much  denser,  has  a  higher  collision 
frequency  and  contains  much  larger  gradients  in  electron  density.  No 
mxnerlcal  results  seem  to  be  available  for  this  type  of  plasma. 

The  basic  references  used  for  the  fundamental  theory  of  plasma 
wave  propagation  were  the  works  by  Bohm  and  Gross  (Ref.  3),  Oster 
(Ref.  4),  Stlx  (Ref.  7),  and  Frledlander  (Ref.  8). 

Three  principal  papers  give  experimental  data  verifying  the 
existence  of  a  propagating  acoustic  wave  in  a  homogeneous  plasma:  Van 
Hoven  (Ref.  9),  Derfler  and  Simonen  (Ref.  10),  and  Malmberg  and  Wharton 
(Ref.  11).  These  papers  basically  verify  the  dispersion  relation  for 
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the  propagating  waves  and  give  experimental  evidence  of  Landau  damping 
of  the  waves. 

Much  work  has  been  done  on  theon.es  for  wave  propagation  In  com¬ 
pressible  plasmas.  Wait  (Ref.  12)  discusses  the  radiation  from  sources 
immersed  in  a  compressible  plasma.  He  shows  that  the  compressibility 
effects  of  the  plasma  on  the  radiation  of  electromagnetic  waves  are 
increased  as  the  ratio  of  the  acoustic  velocity  to  the  speed  of  light 
is  lncressed.  Wait  (Ref.  13)  discusses  the  Influence  of  bouidarles 
oa  wave  propagation  In  a  compressible  plasma.  He  also  describes  the 
absorptive  (impedance)  boundary  condition  for  acoustic  waves.  Felsen 
(Ref.  14)  derives  first  order  coupled  wave  equations  for  the  quasl- 
electromagnetlc  and  quasl-dynamical  fields  unich  he  describes  in  the 
paper.  Burman  (Ref.  15)  derives  alternative  first  order  coupled  wave 
equations  for  the  full-wave  variable  field  quantities  (see  Chapter  III) . 
In  another  paper  (Ref.  5)  Burman  presents  a  detailed  derivation 
of  the  first  order  coupled  equations;  He  also  discusses  the 
coupling  and  power  flow,  and  proposes  some  approximate  solutions. 

Burman  (Ref.  16)  derives  the  second  order  coupled  wave  equations  and 
specializes  the  equations  for  a  planar  geometry.  The  fields  in  a 
region  of  coupling  are  investigated.  Several  approximate  techniques 
for  treating  the  equations  are  discussed.  The  approximate  techniques 
are  not  valid  for  the  type  of  plasma  considered  here. 

When  considering  wave  propagation  in  a  bounded  plasma,  the 
boundary  conditions  become  a  fundamental  problem.  Several  authors  have 
treated  the  subject  of  boundary  conditions  for  a  compressible  plasma. 

Wait  (Ref.  17)  demonstrates  that  the  boundary  condition  which  requires 


die  vanishing  of  die  normal  component  of  the  electron  velocity  leads 
to  results  which  are  consistent  with  cold  plasma  theory.  Yeh  (Ref.  18) 
and  Wait  (Ref.  19)  discuss  the  boundary  conditions  at  a  dielectric- 
plasma  interface,  Sancer  (Ref.  20)  discusses  the  boundary  conditions 
required  for  a  unique  solution  to  the  equations  of  the  hydrod|ynaaic 
treatment.  The  conversion  of  electromagnetic  waves  into  longitudinal 
plasma  waves  at  a  dielectric-plasma  boundary  is  considered  by  Uessel, 
Marcuvitz  and  Shmoys  (Ref.  21). 

A  number  of  authors  have  treated  the  effects  of  the  compressi¬ 
bility  of  a  plasma  surrounding  an  antenna.  Wait  (Ref.  22)  states 
that  it  is  possible  for  a  considerable  portion  of  die  power  to  be 
radiated  as  an  electron  acoustic  type  wave.  The  distribution  of  the 
power  depends  on  the  boundary  conditions.  Kuehl  (Ref.  23)  calculates 
the  radiation  resistance  of  a  short  antenna  in  a  warm  plasma.  Re  uses 
both  a  kinetic  treatment  and  a  hydrodynamic  treatment  and  compares  the 
two  methods. 

A  brief  siamaxy  of  the  topics  covered  in  this  thesis  includes 
Chapter  II  which  gives  a  simplified  derivation  of  the  dispersion 
relation  for  a  pressure  wave  propagating  in  a  compressible  plasma. 

The  plasma  model  is  described  in  detail  and  the  basic  equations  are 
derived.  Chapter  XII  derives  the  second  order  coupled  wave  equations. 
Consideration  is  given  to  both  vertical  and  horizontal  polarizations 
of  Incident  electromagnetic  waves.  In  Chapter  IV  finite  difference 
mmerical  solutions  are  presented  for  the  coupled  wave  equations.  A 
linearly  increasing  plasma  profile  (electron  density  and  temperature) 
with  a  constant  collision  frequency  is  considered.  A  linearly 
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increasing  electron  density  which  increases  to  a  peak  velue,  then  drops 
abruptly  to  zero  is  one  model  considered.  The  ether  model  is  a  linear 
increase  to  a  peak  electron  density  with  a  transition  at  this  point  to 
a  semi-infinite  plaama  which  is  homogeneous  with  an  electron  density 
equsl  to  the  peak  plasma  electron  density.  The  homogeneous 
region  is  considered  semi-infinite  because  of  the  extremely  small 
depth  to  which  acoustic  waves  would  penetrate  this  region  (assuming  a 
colllaion-douineted  plasma) .  Actually  the  homogeneous  region  would  be 
of  finite  extent  cad  at  some  distance  beyond  the  transition  point  the 
electron  density  would  dsersase  to  zero,  i.e.,  a  plasma- free  space 
boundary.  Convergence  of  numerical  solutions  is  consldsred  and 
numerical  solutions  are  compared  with  approximate  analytical  solutions. 
In  Chapter  V  the  eonvmrsl'vt  of  an  electromagnetic  wave  into  e  plasma 
wave  (and  vies  versa)  at  a  dielectric-plasma  intarface  la  studied  as 
s  function  of  the  plasma  properties  and  the  par^ietere  of  the  incident 
electromagnetic  wave. 


CHAPTER  II 


PLASMA  MODEL  DESCRIPTION 

Plasma  Model  Geometry  and  Parameters 

The  planar  geometry  of  either  plasma  model  it,  illustrated  In  Fig. 
1  which  also  pictures  the  Cartesian  coordinate  system  used  throughout. 
The  boundaries  axe  planes  of  infinite  extent  in  the  y  and  z  coordinate 
directions.  The  plasma  parameters  are  as  a  used  to  vary  in  the  x  direc¬ 
tion  only.  The  plasma  can  then  be  described  as  a  layer  (inhomogeneous 
in  the  x  direction)  in  physical  contact  with  a  dielectric  half-space. 
Plasma  parameters  of  primary  interest  include:  electron  density^,  the 
electron-neutral  particle  collision  frequency,  electron  temperature 
and  the  plasma  layer  thickness .  Electron  density  and  temperature  vary 
with  distance  x.  The  collision  frequency  is  assumed  to  be  a  constant 
throughout  the  layer.  Plasmas  of  interest,  in  general  terms,  have 
very  large  peak  (maxlmw)  electron  densities,  very  large  collision 
frequencies  (i.e.,  collision-dominated  plasmas),  and  temperatures 
consistent  with  "warm"  plasmas  (approximately  4000°K). 

Basic  Assumptions 

First  consider  the  regions  surrounding  the  plasma  layer.  The 
dielectric  ic  uniform  and  lossless  with  a  relative  dielectric  constant 
K^.  The  region  outside  the  plasma  layer  is  free  space.  The  permea¬ 
bility  is  that  of  free  space  in  all  three  regions. 

A  plasma  can  be  defined  as  a  gas  containing  a  certain  density 
of  free  positive  and  negative  charges.  It  is  a  well  established 


9 


11 


fact  that  oa  electromagnetic  wave  will  interact  with  charged  particles. 
For  the  electromagnetic  wave  frequencies  and  plasma  parameters  con¬ 
sidered,  the  interaction  produces  organized,  steady -state  oscillations 
in  the  plasma.  Consider  the  frequency  to  be  high  enough  so  that  the 
lens  remain  stationary  and  only  the  electrons  Interact  with  the 
electromagnetic  wave.  Also,  model  the  plasma  as  a  perfect  electron 
gas  with  a  uniform  background  charge  of  positive  ions  making  It 
macro* copl cally  neutral.  It  Is  a  good  approximation  to  consider  die 
ions  os  s  uniformly  smeared  background  charge  for  wavelengths  ouch 
greater  than  the  icteriooic  spacing  (Kef.  3).  Since  this  is  also  a 
requirement  for  any  physical  organized  oscillation  to  exist,  this 
will  always  be  true.  The  assumption  is  made  that  the  wavelengths 
considered  must  be  much  greater  than  the  interionic  spacing  (n^1/3, 
where  represents  the  equilibrium  electron  density). 

The  amplitudes  of  the  fields  in  the  plasma  are  assuned  to  be 
small  enough  so  that  the  equations  can  be  linearized  (Ref.  11) .  This 
means  that  the  eleettic  field  strength,  magnetic  field  strength, 
ordered  electron  velocity,  pressure  field  strength,  and  perturbed 
electron  density  are  considered  first  order  perturbation  quantities. 

Use  an  equilibrium  (zero  order)  electron  density  of  the  plasma 
which  Is  Inhomogeneous  only  In  the  x  direction.  The  plasma  is  con¬ 
sidered  to  be  compressible  and  therefore  has  a  finite  temperature. 
Itself  a  function  of  x.  The  electron-neutral  particle  collision 
frequency  is  assuned  to  be  energy-independent,  independent  of  the 
ordered  electron  velocity,  and  constant  across  the  plasma  layer. 


Hie  plasma  ie  collision-dominated,  l.e.,  v  of  tfae  order  of  u.  This  1« 
e  good  asiiaptia'  for  reentry  type  plasmas. 

The  plasma  is  assimd  to  be  isotropic  vlth  no  external  static  mag¬ 
netic  or  electric  fields.  The  plasma  is  asstmmd  to  be  stationary, 
l.e. ,  the  drifting  velocity  is  aero.  Mo  external  static  body  forces 
actirg  on  the  particles  in  the  plasea  ere  considered. 

The  two  transport  phenomena  of  viscosity  end  heat  conduction  are 
neglected  in  the  plasma.  The  approximation  then  made  la  to  neglect 
the  tens  in  the  equations  of  notion  due  to  viscosity  end  beat  conduc¬ 
tion  (Ref.  8).  The  absolute  error  In  the  approximate  solutions  la  of 
the  order  of  the  product  of  the  distance  through  which  a  wave  propa¬ 
gates  and  of  the  ratio  of  the  tana  neglected  in  the  equations  of 
notion  to  those  retained  (Raf .  8).  Since  the  distances  considered  for 
propagating  waves  in  this  case  an  very  snail,  the  approximation 
should  be  reasonable. 

It  Is  assumed  that  the  plasma  Is  sufficiently  dense  so  that  a 
description  In  tarns  of  a  fores  due  to  e  gradient  of  the  perturbed 
pressure  field  Is  valid.  This  approach  Is  taken  in  solving  plasma 
problems  In  the  ionosphere  (Ref.  24).  The  plasma  presently  considered 
Is  much  denser  than  the  ionospheric  plasma,  and  this  assunptlon  thus 
a  more  valid  one. 

There  are  two  important  phenomena  in  acoustic  wave  propagation 
which  do  not  appear  in  the  hydrodynamic  treatment.  Physically,  a 
plasma  oscillation  with  a  wavelength  shorter  than  the  inter-particle 
spacing  is  meaningless  in  terms  of  collective  motion.  This  has  been 
discussed  earlier.  The  other  phenomenon  is  a  r.  on  collision  al  damping 
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effect  celled  Lender  deplo|  (Ref.  25).  The  organized  plan*  os  ci  na¬ 
tions  ate  damped  by  e  treiefci  ot  energy  free  the  wave  to  the  energy 
la  the  randan  thermal  motions  cf  electrons  which  are  moving  et  e  speed 
dose  to  the  wave  velocity.  Landau  damping  becomes  important  when  the 
acoustic  wavelength  approaches  the  Debye  length  of  the  plasma.  When 
the  acoustic  wavelength  is  smaller  than  die  Debye  length,  collective 
motion  ie  essentially  destroyed  (Ref.  3).  It  will  be  ass  used  that  in 
regions  where  a  pressure  field  (acoustic  wave)  exists,  the  acoustic 
wavelength  will  be  larger  then  die  Debye  length;  however,  since  the 
plasma  is  inhomogeneous ,  an  acoustic  wavelength  is  difficult  to  deter- 
mine.  The  asauaption  will  be  verified  in  Appendix  D  by  calculating 
the  wavelength  from  the  dispersion  relation  for  a  lossy,  homogeneous 
plasma  having  an  electron  density  equal  to  the  electron  density  at 
discrete  points  along  the  inhomogeneous  plasma  profile. 

It  is  necessary  now  to  specify  the  frequency  range  of  interest. 

The  constraint  on  the  plasma  wavelength  due  to  Landau  damping  deter¬ 
mines  a  maximal  frequency  of  operation.  For  a  collisionless,  homoge- 
eous  plasma,  Bohn  and  Gross  (Ref.  3)  have  determined  the  maximum 
frequency  to  be  approximately  *2  fp,  where  f  is  the  plasma  frequency. 
Rather  than  specifying  a  maximum  frequency,  the  approach  in  the  pre¬ 
vious  paragraph  will  be  used  and  frequencies  larger  than  the  frequency 
determining  an  acoustic  wavelength  equal  to  the  Debye  length  will  not 
be  considered. 

The  electron  plasma  frequency  Is  basically  a  cutoff  frequency. 

Below  the  plasma  frequency,  the  electrons  can  respond  to  an  electro¬ 
magnetic  field  but  above  it  they  cannot  follow  the  field.  Actually, 
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Che  cutoff  is  not  that  sharp.  Van  Howes  (Ref.  *)  states  that  below 
the  plasma  frequency,  acoustic  waves  are  strungiy  dasped  by  collisiosal 
effects.  It  is  assumed  that  the  hydrodynanic  approach  (including 
collisonal  damping)  can  be  used  to  investigate  aCO'wtic  wave  phenomena 
at  frequencies  well  below  the  plasma  frequency.  There  is  no  reason 
why  the  electron  collective  notion  due  to  the  interaction  between  the 
electromagnetic  field  and  the  charged  particles  should  not  stll^  be 
described  by  the  hydrodynaadc  equations,  since  die  electrons  can  still 
collectively  interact  with  the  electromagnetic  field  for  frequencies 
below  die  electron  plasma  frequency.  A  minimum  operating  frequency 
would  be  the  ion  plasma  frequency,  at  which  ion  motion  would  have  to 
be  considered. 

Rationalized  MRS  units  axe  used  for  all  n tarn rl cal  calculations. 

Basic  Equations 

The  basic  equations  are  the  two  Maxwell's  curl  equations  and  two 
equations  from  hydrodynasdc  theory  describing  the  dynamics  of  a  con¬ 
tinuum.  These  equations  fora  the  aathenatl cal  model  for  calculating 
the  electromagnetic  and  acoustic  fields  in  die  plasma. 

Maxwell's  equations  are  given  below,  where  it  is  assumed  that  the 
regions  of  space  of  interest  are  located  far  enough  away  from  any  mag¬ 
netic  or  electric  sources  so  that  their  effect  is  negligible. 

7  x  i  If  <» 

a?  -> 

•  c  +  (n  +  n)  eV 
o  9t  o 


v  x  3 


(2) 
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A  convection  current  density  ten*  has  been  included  in  the  second 
equation.  The  equillbriias  electron  density  is  given  by  and  the 
electron  density  perturbation  (fluctuation)  is  given  by  n.  The  charge 
accumulation  or  depletion  is  given  by  ne  which  is  in  general  a  function 
of  both  time  and  position.  The  dbarge  of  a  single  electron,  e,  is  a 
negative  nimber.  The  electron  density  fluctuation,  n,  can  be  either 
positive  or  negative.  The  acoustic  vector  field,  tf,  is  the  ordered 
velocity  of  the  electrons. 

The  equation  for  conservation  of  momentum  for  an  in  viscid,  non- 
heat-conducting  acdU.ua  is  given  by  Frledlander  (Ref.  8).  This  equation 
is  soaetlaea  called  the  hydrodynaedc  equation  or  dynastic  equation. 

(»o  *“)  | It  +  5}+  v(po ♦ »)  ‘  •(”<> +  “) 1 

♦  «(?  X  »*)  („o  +  n) 

-  («„  +  <•)  •  <3> 

The  equilibrium  quantities  are  si&scripted  and  the  perturbation 
quantities  are  not.  The  ness  density  is  given  by  p  and  is  equtl  to 
the  mass  of  an  electron  times  the  electron  number  density.  The  colli¬ 
sion  frequency  is  given  by  v  and  its  effect  is  included  in  the  momentum 
equation  as  a  drag  force  proportional  to  the  ordered  velocity  of  the 
electrons. 

Equation  (3)  contains  only  one  zero  order  term,  VpQ.  Friedlander 
(Ref.  8)  has  shown  that  for  a  plasma  at  rest  with  no  static  body  forces 
acting  upon  it,  VpQ  •  0  even  if  the  plasma  is  inhomogeneous.  This  is  a 
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consequence  of  hydrostatic  eqid’roriiflo.  Since  pQ  «  n  mV*  this  con¬ 
dition  would  imply 

%  •  '("oS1)  •  ST™o  +  SV1  ■  0  •  “> 

where  Kg  i»  Bultzmcnr.’s  constant  and  T  is  the  plasma  temperature.  This 
imposes  a  definite  relationship  between  nQ  and  T.  It  is  not,  however, 
the  relationship  which  has  been  assumed,  i.e.,  that  the  temperature 
profile  follows  the  electron  density  profile  in  form.  This  type  of 
temperature  variation  is  assumed  because  it  is  a  good  approximation 
to  the  temperature  profile  in  a  reentry  type  plasma  where  higher 
temperatures  result  in  more  ionization  and  give  a  higher  value  for 
nQ.  It  is  concluded  that  VpQ  4  0.  This  case  is  discussed  by  Burman 
(Ref.  16).  Equation  (3)  must  now  be  modified  to  Include  a  zero  order 
term  which  will  cancel  the  VpQ  term,  since  VpQ  4  0.  The  modification 
is  written: 

(Po  +  c)  (f!+ 1?-”  *}+  7(p0  + »)  -  (»o  +  “)  (K +  f) 

■  *(”o  +  ")  1  +  e(?  *  W0S)  (n0  +  n) 

-  (»o  +  “)  ^  •  «*> 

where  p  ?  is  the  zero  order  term  required  to  balance  the  zero  order 
o  o 

pressure  gradient.  A  first  order  term,  p  f  +  ,  is  also  introduced. 

o  o 

The  force  term  (p  +  p)(?  +  ?)  is  not  due  to  any  external  static  fields 
oo 

or  forces  since  these  are  assumed  negligible,  but  is  due  to  an  ambi- 
polar  diffusion  electric  field  due  to  the  static  pressure  gradient. 
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The  zero  order  part  of  this  force  tern  balances  the  static  pressure 
gradient  while  the  effect  of  the  first  order  part  is  considered  to  be 
negligible  in  the  first  order  part  of  the  momenta*  equation.  The  first 
order  conservation  of  moaentus  equation  then  becomes : 

■(”„  +  »)  | !?  +  <*•’>  V  j  +  9p  -  «(n0  +  r)l  +  «(v  X  „o«)  (no  +  a) 

+  nj  vV  .  (6) 

The  exact  equation  of  conservation  of  mass  Is  given  by  the  con¬ 
tinuity  equation 

M°o  +  ‘,)+'-{(‘>o  +  C)  (?o  +  ’)|‘°  •  (7) 

The  vector  Vq  represents  any  drifting  velocity  of  the  plasma  and  has 
been  assumed  to  be  zero.  In  addition,  since  the  equilibrium  mass 
density  is  not  a  function  of  time,  Eq.  (7)  can  be  simplified  as  follows: 


It  is  convenient  at  this  point  to  linearize  Eq.  (8).  This  is  done  by 
neglecting  the  second  order  term  pV  inside  the  brackets.  As  stated  in 
the  basic  assumption,  the  amplitudes  of  the  disturbances  are  small 
enough  so  that  the  linearization  approximation  is  valid.  Then 

If  +  v.  p0V  -0  .  (9) 

Applying  the  vector  identity  ?•(<$)  *  <p  +  £*V<J>,  one  can  write 
Eq.  (9)  as 
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m  ~  «  7*V  +  mV* Vn  -  0  .  (10) 

at  O  O 

*♦ 

The  term  V-v’p^  is  neglected  in  the  majority  of  works  on  acoustic  wave 
propagation  (Refs.  5,  14,  IS,  and  16).  This  is  valid  for  a  homogeneous 
plasma  (since  7pQ  ■  0);  however,  this  asswption  may  not  be  valid  for 
a  plasma  with  steep  electron  density  gradients.  For  generality  this 
term  will  be  retained  in  this  derivation  of  the  wave  equations. 

Now  combine  Eq.  (10)  with  the  approximate  adiabatic  equation  of 
state  for  a  perfect  electron  gas: 

<n) 

where  y  is  the  ratio  of  specific  heat  at  constant  pressure  to  specific 
heat  at  constant  voluse.  The  quantity  y  is  normally  ass  used  to  be  3 
for  an  electron  gas  (Ref.  25).  For  a  perfect  electron  gas,  take 

Po  -  no!^T  (12) 

and 

P  «  nKgT  (13) 

where  Kg  is  Boltzmann's  constant  and  T  is  the  temperature. 

Multiplying  both  sides  of  Eq.  (10)  by  y/no  and  using  Eqs.  (11), 
(12) ,  and  (13)  yields 


lE 

at 


+  mn 


-  /yKb  T\. 

7*V  +  m  v‘VnQ 


0  . 


(14) 


The  term  in  parenthesis  is  the  square  of  the  acoustic  velocity  (rms 
thermal  velocity)  in  the  electron  gas. 
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The  final  font  of  the  equation  ia  then 

I®*  +  m  u^V*^  +  at^V^Tn  ■  0  . 

dt  o  o  o  o 


(16) 


How  aasiae  that  all  of  the  field  quantities  have  a  steady-state 
time  dependence  expressed  by  the  factor  e^ut.  This  factor  will  be 
suppressed  and  the  field  quantity  aaplltudes  will  be  considered  as 
phasors  having  coop  lex  values  which  are  f  Actions  of  position  only. 
Then  If  3/3t  Is  replaced  by  jut,  Eqs.  (1) ,  (2),  and  (6)  are  linearized, 
and  Eq.  (16)  is  used: 


V  X  2  -  ,  (17) 

7XB*  juieQ  t  +  nQ  eV  ,  (18) 

an  (Ju  +  v)$  *  n  et  -  Vp  ,  (19) 

o  o 

jwp  +  an  u2V*V  +  mv^V'Vn  *  0  .  (20) 

O  O  0  0 


Compressibility  and  the  Propagation  of  Pressure  Waves 


It  will  now  be  shown  that  an  acoustic  or  pressure  field  can 
exhibit  wave  behavior  only  In  a  compressible  plasma.  For  simplicity, 
assune  a  plasma  that  is  homogeneous,  lossless  and  of  Infinite  extent. 

The  plasma  Is  excited  by  a  steady-state  oscillator  at  some  point  which 
Is  far  from  the  observation  region.  Equilibrium  quantities  are  sub¬ 
scripted  by  a  zero.  The  first  order  perturbations  are  not.  The  basic 
equations  were  taken  from  Stlx  (Ref.  7)  and  specialized  to  a  one- 
component  electron  plasma.  The  equations  are  linearized  by  neglecting 
products  of  all  first  order  terms  since  small  perturbations  are  assumed. 
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First  consider  on  inexpressible,  or  cold, 
temperature. 


po  “  W  “  0 

and 


p  »  nKgT  »  0 


The  con  wet  ion  current  density  i»  given  by 


3  “  n  eV  . 
o 


The  conservation  of  moment!*  equation  is 


3V  ± 

m  *  e£  , 


since  p  »  0. 

The  conservation  of  charge  equation  is 


ii  + 

3t 


where  the  charge  density  is  p  ■  en. 
From  Gauss'  law. 


■  p/e  ■  en/e 
o  o 

Taking  the  partial  derivative  of  Eq.  (25)  with 

at2  3t 


plasma  with  aero 

(21) 

(22) 

(23) 

(24) 

(25) 

(26) 

resect  to  time  gives 
(27) 


Taking  the  partial  derivative  of  Eq.  (23)  with  respect  to  time  and 
using  Eq.  (24)  gives 
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-»  e£n 

3J  m  o 
3t  m 


(28) 


Sihstitutlng  Eq.  (23)  into  Eq.  (27)  and  using  Eq.  (26)  gives 

32n 


e2n 


at 


o 

E 

O 


+ - -  n  -  0 

2  ■*- 


or 


--~n-  +  «2  n  ■  0  , 

at2  ? 


(29) 


where 


(30) 

V  “o 

is  the  electron  plasma  frequency. 

Although  the  spatial  dependence  of  a  pressure  field  cannot  be 
determined  in  this  case  (since  p  *  0),  acoustic  phenomena  can  be 
related  to  charge  fluctuations,  since  the  acoustic  wave  results  from 
the  collective  motion  of  the  particles.  From  Eq.  (29)  it  is  observed 
that  the  collective  motion  of  the  electron  density  perturbation  (or 
charge  fluctuation  since  p  -  en)  is  oscillatory  in  time  but  is  spa¬ 
tially  independent.  This  disturbance  does  not  propagate. 

For  a  compressible  plasma,  Eqs.  (23),  (25),  and  (26)  still  apply 
but  a  pressure  term  must  be  included  in  the  conservation  of  momentum 
equation . 


$  ■* 
mn  r—  «  en  E  -  Vp 
o  3t  o  r 


(31) 


For  a  finite  temperature, 
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’o  *  noV 

and 


p  -  nKjT  . 


The  approximate  adiabatic  equation  of  atate  for  an  ideal  gas  is  given 
by 


Using  Eqs.  (23),  (23),  (26),  and  (31)  and  the  aiDe  procedure  used 
previously  gives 


■  +  u2n  -  —  ?2n  ■  0 
3t2  ?  ■ 


or 


u2  V2p  -  0  , 


(32) 


where  u^  is  the  acoustic  velocity  (rms  thermal  velocity)  In  the  elec¬ 
tron  gas  and  is  given  by  Eq.  (15). 

For  steady-state  tine  dependence  e^*,  p  has  the  form 
+  P2*  (^Wt+^P*r),  where  pj  and  p2  are  constants.  This  is  the  mathe¬ 
matical  representation  for  a  propagating  plane  wave.  If  this  form  of 
p  Is  substituted  into  Eq.  (32),  the  dispersion  relation  is  found  to  be 


u)2  ■  u)2  +  u2  K2  .  (33) 

pop 

It  is  concluded  that  a  pressure  wave  can  propagate  in  a  compressible 
plasma  with  a  wavelength  (X^  ■  2it/K^)  determined  by  Eq.  (33)  for  & 


given  frequency.  This  is  Identical  to  the  equation  obtained  by  Bohm 
and  Grass  (Ref.  3). 

For  ao  inhomogeneous  plasma,  it  is  not  possible  to  describe  the 
acoustic  vase  in  tares  of  a  unique  propagation  constant.  It  is  also 
fouid  that  an  acoustic  wave  is  coupled  with  an  electromagnetic  wave  in 
an  inhomogeneous  plasma.  The  propagation  of  acoustic  waves  in  an 
inhomogeneous  plasma  is  described  by  the  coupled  wave  equations  derived 
in  the  next  chapter. 


CHAPTER  III 


DERIVATION  01  THE  COUPLED  WAVE  EQUATIONS 
Full-Wave  Variable* 

In  a  homogeneous  plasma,  the  electromagnetic  field*  and  die  dynamic 
(plasma)  field*  can  be  aeparated  into  two  distinct  modes.  Energy 
exchange  between  the  two  modes  can  occur  only  at  an  interface  and  each 
mode  propagates  independently  of  the  other.  There  then  exists  a  pure 
electromagnetic  field  and  a  pure  dynamic  (plcsma)  field.  The  electro¬ 
magnetic  field  contains  all  of  the  magnetic  field  and  the  dynastic  field 
contains  all  of  the  charge  accimulatlon  (Ref.  14). 

In  an  lnhomogensous  plasma  (Ref.  6)  the  separation  of  the  various 
fields  into  pure  electromagnetic  or  pure  dynamic  fields  is  impossible 
due  to  the  spatial  variation  of  the  plasma  properties .  Thus  an  electro¬ 
magnetic  field  quantity  would  consist  of  an  electromagnetic  component 
plus  a  dynamic  component  and  therefore  could  not  be  truly  a  transverse 
wave.  Similarly,  a  dynamic  field  quantity  would  consist  of  a  dynamic 
component  plus  an  electromagnetic  component  and  therefore  could  not  be 
truly  a  longitudinal  wave.  This  is  a  fuidamental  description  of  the 
coupling  between  the  two  field  quantities.  It  is  difficult,  however, 
to  investigate  quantitatively  the  coupling  effects  by  evaluating  the 
electromagnetic  component  and  the  dynamic  component  in  a  given  field 
quantity.  These  components  are  referred  to  as  a  quasi-electromagnetic 
mode  and  a  quasi-plasma  mode  (Ref.  14).  Since  the  given  field  quanti¬ 
ties  are  described  physically  by  the  sun  of  these  quasi-modes,  it  is 
desirable  to  obtain  solutions  in  terms  of  the  total  field  quantity. 
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i.e.,  the  sub  of  the  Bodes.  These  ere  referred  to  in  the  literature 
st  full-wave  variables  (Kef.  6) ,  and  are  simply  the  field  quantities 
described  by  die  basic  equations  given  in  Eqs.  (17),  (18),  (19),  and 
(20).  The  wave  equations  describing  wave  propagation  in  an  iuhonoge- 
neous  plasaa  are  derived  for  the  full-wave  variables  using  the  two- 
diaenalonal  geometry  of  Fig.  1. 

Vertical  Polarisation 

The  electromagnetic  and  pressure  fields  are  shown  in  Fig.  2  for 
a  vertically  polarized  Incident  electromagnetic  wave.  The  electro¬ 
magnetic  field  quantities  £  ana  £  are  vectors.  The  pressure  p  is  a 
scalar  quantity  represented  In  Fig.  2  as  a  wave  propagating  at  some 
angle  0^.  A  steady-state  situation  Is  assused  with  all  field  quan¬ 
tities  having  an  e^ut  time  dependence.  The  radian  frequency  of  the 
electromagnetic  wave  is  u.  The  *^ut  factor  will  be  suppressed  in 
die  derivation  since  It  is  common  to  all  field  quantities. 

Since  the  plasaa  Is  homogeneous  In  the  z  direction  and  the 
Incident  wave  Is  a  plane  wave,  there  is  no  loss  of  generality  in 
assuming  that  the  fields  do  not  vary  in  the  z  direction  and  3/3z  »  0 
for  any  field  quantity. 

The  electric  field  has  components  E_  and  E^.  The  magnetic  field 

has  only  a  z  component,  H  .  The  velocity  (ordered  electron  velocity) 

field  has  components  V  and  V .  The  pressure  field  has  only  a  magni- 

*  / 

tude  p.  These  are  all  complex  phasor  quantities  containing  both 
amplitude  and  phase  information.  They  are  also  fmctions  of  the 
spatial  coordinates  x  and  y.  Instantaneous  values  are  found  by 
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multiplying  by  e^wt  and  taking  the  real  part. 

the  linearised  basic  equations  are  given  again  for  reference. 


7XE-  -jwuo  H  (34) 

7  2  H  «  jut  E  +  neV  (35) 

o  o 

an  (Jji)  +  v)V-neE-Vp  (36) 

o  o 

j«p  +  «&  u £  V»?  +  au2  V*7n  -  0  (37) 

o  o  o  o 


Equations  (34)  to  (37)  will  now  be  written  in  component  form  using 
the  field  coop oaents  listed  above,  the  approach  used  is  to  solve  for 
all  of  the  field  components  in  tens  of  H  and  p.  The  wave  equations 
will  be  derived  for  and  p  and  the  solutions  will  allow  all  of  the 
fields  in  the  plasa*  to  be  determined. 

Froa  Eq.  (34), 


3E  3E 

— — s—  —  -liuU  H 
3*  sty  J  s 


(38) 


Froa  Eq.  (35), 


3E 

— r—  -  Jue  E  +  n  e  V  , 
3y  J  o  x  ox* 


(39) 


and 


From  Eq.  (36), 


3H 
_ z 

3x 


j^0 


E 

y 


+  n  e  V 
o  y 


(40) 


mn  (jo)  +  v)  V 
o  J  x 


_  l£ 


en  E  . 
ox  3x 


> 


(41) 
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am  (jw  +  v)  V  -  en  E  -  . 

o  y  o  y  3y 


From  Eq.  (37), 


(av  av  \ 

r~~  +  ) 

3  x  3y  / 


3n 

+■  nu2  V  — — 
o  x  3z 


For  simplicity,  define 


A  ■  ju  +  v  , 


U  -  1  -  J  v/w  . 


Since  the  collision  frequency  has  been  assisted  constant,  neither  A  nor 
U  are  functions  of  the  spatial  coordinates. 

Solving  Eq.  (41)  for  and  si&stltutlng  this  Into  Eq.  (39)  gives 


Jaj£  fi.i)B  +  . 

0  \  w2U/  X  3y  Am  3x  * 


where  u>2  has  been  defined  in  Eq.  (30).  Now  define  the  two  quantities 


e  ■  e  e 

o  p 


w 

i-4-  . 


where  c  Is  the  complex  permittivity  of  the  plasma. 
Equation  (44)  can  then  be  solved  for  E^  to  give 


t  on  . 

E  -  -i - 1+  -.e-  32. 

x  juie  3y  jweAo  3x 


Solving  Eq.  (42)  for  and  substituting  this  into  Eq.  (40)  gives 

E  .  _  _i _ l  +  e  l£ 

y  joje  3x  ju>eAm  3y 


(48) 
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Solving  Eq.  (41)  for 


Solving  Eq.  (42)  for  Ey 


and  substituting  this  into  Eq. 

c*H 


_  O 

x  4b>cAa  3y  cAmn^  °* 


V  -- 


and  substituting  this  into  Eq. 


au  c 

«■•  __£  _  0  j>£ 

juitAa  ~3x  *  cAaa^  3y 


(19)  gives 

(49) 

(40)  gives 

(50) 


Equations  (47)  to  (00)  agree  with  those  obtained  by  Walt  (Ref.  22). 

It  Is  concluded  that  all  field  quantities  in  the  plasma  can  be  deter¬ 
mined  if  aolutlons  for  H  and  p  are  found. 

z  v 

Since  the  plasma  is  homogeneous  in  the  y  direction,  the  y  compo¬ 
nent  of  Che  propagation  constant,  denoted  by  K^,  is  conserved  for  all 
x  values.  The  y  component  of  the  propagation  constant,  K^,  of  the 
plane  electromagnetic  wave  in  the  dielectric  is  given  by 


(51) 


where  Kd  is  the  dielectric  constant.  Following  references  5,  12,  13, 
15,  and  16,  the  assumed  forms  of  the  magnetic  field  and  the  pressure 
in  the  plasma  are 


H 


H2(x) 


j*  y 
,  y 


(52) 


and 


p(x) 


jK  y 

,  y* 


(53) 


All  other  field  components  are  assumed  to  vary  in  the  same  manner. 

The  problem  then  becomes  essentially  one-dimensional,  since  from 
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Eqs.  { 52 >  and  (53), 


K 


z 


and 


i£ 

3y 


P 


describe  the  y  variat’oas  of  H  and  p.  Therefore,  3/3y  can  be  replaced 

z 

by  JXy  and  the  derivatives  with  respect  to  x  be cane  total  derivatives. 

The  component  fores  of  the  basic  equations  will  now  be  written 

with  the  y  dependence  of  the  field  conponents  ass  used  above.  Since  the 

jK  y 

e  y  factor  is  common  to  all  the  field  components,  it  will  be  suppressed. 
dE  (x) 

-  jKy  Ex(x)  -  -jwuo  Hs(x)  (54) 

jKy  Hx(x)  -  j»£o  Ex(x)  +  noe  V%(x)  (55) 

dH  (x) 

- *T"  “  jueo  Ey<x)  +  noe  Vy(x)  (56) 

tmQ  A  Vx(x)  -  enQ  Ex(x)  -  (57) 

ono  A  Vy(x)  -  enQ  Ey(x;  -  JKy  p(x)  (58) 

/dV  (*)  \  dn 

juip (x)  +  mnou2  y  ■  -  +  JKy  Vy(x)j  +  mu*  Vx(x)  -  0  (59) 

The  field  components  are  now  expressed  as  functions  of  x  only.  The 
functional  dependence  on  x  will  be  assumed  and  in  the  following  deri¬ 
vation  the  notation  Hg(x)  will  be  replaced  by  Hz  and  similarly  for  the 
jther  field  components.  Equations  (47)  to  (50)  will  be  used  together 
with  Eqs.  (54)  to  (53)  to  derive  the  coupled  wave  equations. 
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Taking  the  derivative  of  Eq.  (56)  with  respect  to  x  gives 


d*H  dE 

z  .  y  .  o 

- *  JwE  j  +  e  — 

^2  o  dx  dx 


("oV  y) 


where  the  product  nQV^  must  be  retained  in  the  derivative  operation 
since  n  sad  V  are  flections  of  x.  Substituting  dE  /dx  found  from 

o  7  y 

Eq,  (54)  and  n^V^  foiad  from  Eq.  (50)  into  Eq.  (60)  gives 


- -  -  juc  -  juu  H  +  JK  E 

dx2  0  [  o  z  y 


.  r  en  dH  jK  e  1 
d  _  o  _ z  _  J  y  o  I 

dx  ~  JideAjb  dx  ~  eAm  PJ  * 


Substituting  fotnd  from  Eq.  (47)  into  the  equation  above  gives 


d2H 

- -  -  K2H  -  ue  1 

dx2  0  *  ° 


K  [-*  H 

i  y  I  ide 


+  e  dP 

z  jiueAn  dx 


it- 


en  dH  jK  e 

O _ Z  y 

JweAm  dx  eAm 


L°p] 

to  PJ  ’ 


where  K2  -  u2woeo.  Expanding  the  derivative  in  the  third  term  on  the. 
right-hand  side  of  the  previous  equation,  noting  that  nQ  and  e  are 
functions  of  x,  and  performing  the  necessary  multiplications  gives 

d2H  K2  K  e 

- E.K2B  -  -Z  H  -t-V* 

.2  o  z  e  z  je  Am  dx 
dxz  p  J  p 


d2“z .  i  d2"z .  r« dH 


+  Jsr 


eK  eK 

+ _ Z_  _  1  _JL  _V_£f  p 

jtpAm  dx  J  Am  L  dx  J  p 


After  some  simplification; 


d*H 

— -  +  , 

dx2  P 


’d (xb.  Y]  dH  f  1  eK  fdfVe  ) 

_i_E/  _I  +  K2t  -  K2  H  -  — £  c  ^~.Sl  p  . 
dx  dx  Lop  yl  z  wmU  p  dx  J  v 


Since 


'N, 

dx 


i  dc 

.  i _ E 

,-2  dx  * 


the  final  form  of  the  first  of  the  two  coupled  differential  equations 
i,.i  Hz  and  p  is  given  by: 


d*  H 

,  de 

dH 

r 

eK 

,  de 

2 

1 _ E 

__ £  + 

|K2c 

[°P 

-*] 

H _ * 

z  umU 

1 _ E 

dx- 

c  dx 

P 

dx 

c  dx 

P 

The  coupling,  in  this  case,  is  introduced  by  the  term  Involving  p  on 
the  right-hand  side  of  the  above  equation  and  the  coefficient  of  p  is 
called  the  coupling  coefficient. 

The  second  coupled  differential  equation  will  be  obtained  from 
Eq.  (59),  repeated  here  for  reference. 


[dV  1  dn 

-7-2-  +  jK  V  +  mu2  V  -r-£  -  0  . 

dx  y  y  o  x  dx 

Taking  the  derivative  of  (Eq.  (49))  with  respect  to  x,  multiplying 

V  (Eq.  (50))  by  jK  ,  substituting  for  V  and  replacing  these  terms 
y  y  x 

in  the  above  equation  by  the  expressions  just  found  gives. 
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Equation  (62)  is  the  second  coupled  differential  equation  in  Hz  and 
p.  The  coupling,  in  this  case,  is  introduced  by  the  term  involving  H 

z 

on  the  right-hand  side  of  Eq.  (62).  The  two  Eqs.  (61)  and  (62)  are 


Identical  to  those  derived  by  Burman  (Ref.  16)  if  the  collision  fre¬ 


quency  is  assumed  constant  and  the  third  term  on  the  left-hand  side  of 
Eq.  (37)  is  neglected  (which  is  not  necessarily  valid  for  the  plasmas 
considered  here).  When  this  term  is  retained  in  the  derivation,  the 
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0  ■  r.  eV 
o  x 


(66) 
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0  ■  n  eV 
o  y 


(67) 


3H  3H 

-r-2-  -  ~  -  Jwc  E  +  n  eV 
3x  3y  or  ox 


(68) 


3V  3V  3V  dn 

jwp  +  am  u2  I  v~  +  +  -r-~  I  +  mu2  V  -r~  *  0  (69) 

J  K  o  ol  3x  3y  3r  J  o  x  dx 

Since  H  and  H  do  not  vary  in  the  s  direction,  Eqs.  (66)  and  (67) 

y  ^ 

show  that  Vx  and  are  identically  zero  everywhere  in  the  plasma. 
Equation  (68)  indicates  that  there  can  be  a  s  cooponent  of  the  veloc¬ 
ity;  however,  it  has  been  assumed  that  all  field  components  are 

constant  in  the  x  direction  so  3V  /3x  «  0.  Also,  since  V  and  V  are 

*  x  y 

zero  everywhere,  3Vx/3x  •  0  and  3Vy/3y  »  0.  Equation  (69)  then  becomes 

jwp  -  0 

or 

p  ■  0  . 

This  implies  that  the  pressure  field  is  identically  zero  everywhere 
in  the  plasma;  hence,  the  acoustic  mode  is  not.  excited. 


CHAPTER  IV 


NUMERICAL  SOLUTION  OF  THE  COUPLED  WAVE  EQUATIONS 
Coupled  Equations 

The  tern  "coupled  equations'"  describee  a  set  of  simultaneous 
differential  equations  with  the  following  properties  (Ref.  26): 

<l)  There  is  one  independent  variable  which  in  this  derivation 
is  the  distance  x. 

(2)  The  number  of  equations  is  the  same  as  the  number  of  depend¬ 
ent  variables, 

(3)  In  each  equation  one  dependent  variable  appears  in  deriva¬ 
tives  of  higher  order  than  the  others.  The  terms  with  this  variable 
are  called  "principal  terms  and  the  remaining  terns  "coupling"  teems. 

(4)  The  principal  terms  contain  a  different  dependent  variable 
in  each  equation,  so  that  each  dependent  variable  appears  in  the 
principal  terms  of  one  and  only  one  equation. 

Clecmow  and  Heading  (Ref.  27)  add  one  more  property  when  the  equations 
describe  wave  propagation: 

(5)  In  a  homogeneous  mediun  the  right-hand  sides  of  the  equa¬ 
tions  vanish  since  the  coefficients  of  these  coupling  terms  contain, 
as  factors,  derivatives  of  the  properties  of  the  mediun.  The  left- 
hand  sides  then  give  the  characteristic  waves.  Equations  (61)  and 
(62)  have  all  of  the  properties  described  above. 

There  are  several  methods  which  can  be  used  In  solving  coupled 
wave  equations.  Solutions  can  be  found  by  successive  approximations 
where  the  first  approximation  is  to  neglect  the  coupling  terns  and 
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solve  the  resulting  equations.  The  values  obtained  for  the  dependent 
variables  are  substituted  in  the  coupling  terns  and  the  resulting 
equations  are  solved  to  give  a  better  approximation.  Where  a  coupling 
tern  is  large,  the  method  of  successive  approximations  cannot  be  used. 
When  the  coupling  can  be  neglected,  WKB  solutions  can  be  obtained; 
however,  these  solutions  fail  near  regions  of  reflection  or  coupling 
(Ref.  26).  The  two  equations  can  be  tsacoupled  by  combining  the  two 
into  a  single  fourth  order  differential  equation  in  one  dependent 
variable.  This  equation  must  then  be  solved.  The  two  second  order 
equations  cm  also  be  reduced  to  four  first  order  equations  which  are 
also  coupled.  These  can  be  solved  by  matrix  methods.  Burra  an  (Ref.  15) 
gives  four  different  methods  of  solution  for  the  cowled  wave  equations; 
however,  s  slowly  varying  medium  is  asswed.  Another  possible  solution 
is  to  integrate  numerically  the  coupled  wave  aquations.  This  is  the 
method  used  later  in  this  chapter.  The  main  advantage  of  this  method 
Is  that  coupling  terms  ate  not  neglected  and  the  solutions  are  valid 
in  media  with  rapidly  varying  properties. 

The  existence  and  uniqueness  of  solutions  obtained  for  Eqs,  (17) 
to  (20)  are  considered  here.  Sancer  (Ref.  20)  has  discussed  the  bound¬ 
ary  conditions  which  yield  unique  solutions  to  the  linearized  warm 
plasma  equations  which  are  used  in  this  derivation.  Solutions  to  these 
equations  exist,  for  simple  geometries,  if  one  requires  that  2  +  M 
scalar  boundary  conditions  are  satisfied  at  each  boundary.  The  nunber 
of  different  types  of  particles,  M,  considered  in  the  equations  of 
motion  is  1  for  an  electron  gas.  The  boundary  conditions  used  are  the 
two  electromagnetic  boundary  conditions  (continuity  of  tangential  E 


acd  iJ)  ,  plus  a  single  acoustic  boundary  condition  which  depends  on  the 
type  of  boundary  considersd  Therefore,  solutions  to  the  equations 
can  be  found.  The  two  different  acoustic  boundary  conditions  are  the 
vanishing  of  the  normal  component  of  the  electron  velocity  and  the 
continuity  of  the  pressure  across  a  boundary.  Sancer  shows  that  eithe 
of  these  conditions  when  used  with  the  two  electromagnetic  boundary 
conditions  leads  to  a  unique  solution. 

Plasma  Profiles 

The  plasmas  of  interest  have  steep  electron  density  gradients 
near  the  wall  (x  «*  0)  rising  to  a  peak  value  at  a  distance  very  near 
the  wall.  The  electron  density  then  decreases  from  the  peak  as  the 
temperature  becomes  lower  and  ionization  is  less.  Consideration  is 
given  to  very  thin  plasaa  "sheaths"  where  the  region  of  interest  is 
5  mm  or  less  away  frcm  the  wall.  Electron  density,  temperature,  and 
collision  frequency  profiles  are  illustrated  for  a  typical  sheath  in 
Fig.  4,  The  temperature  profile  is  similar  in  form  to  the  electron 
density  profile.  This  is  due  to  the  fact  that  the  number  of  electrons 
depends  on  the  degree  of  ionization  which  iu  turn  depends  on  the 
temperature.  For  a  reentry  type  plasma,  the  high  temperatures  are 
associated  with  aerodynamic  heating  which  occurs  when  a  vehicle 
reenters  the  atmosphere.  The  collision  frequency  profile  is  almost 
constant  across  the  sheath.  As  a  first  approximation,  the  collision 
frequency  will  be  taken  to  be  constant.  The  collision  frequency  will 
be  larger  for  larger  peak  electron  densities. 

An  electron  density  profile  that  has  a  linear  increase  from  the 
wall  to  the  peak  is  used.  This  profile  is  simple  enough  for  a  fairly 
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straight  forward  numerical  analysis;  however,  by  taking  the  lirear 
Increase  near  the  wall  large  enough,  a  meaningful  study  of  the  effect 
of  plasma  inhomogeneity  on  pressure  wave  propagation  in  the  plasma  can 
be  made.  There  are  two  basic  rrofiles  considered  with  the  linearly 
increasing  electron  density.  The  first  profile  is  shown  in  Fig.  5(a). 
In  this  case,  the  electron  density  is  assumed  to  decrease  abruptly  to 
zero  at  the  outer  boundary.  The  outer  bowdary  is  then  just  a  plane 
dividing  the  plasma  and  free  space.  The  second  profile  is  shown  in 
Fig.  5(b).  Here  the  electron  density  beyond  the  peak  is  taken  to  be 
a  constant  equal  to  the  peak  electron  density.  In  this  homogeneous 
region,  the  wave  equations  are  uncoupled.  The  electromagnetic  and 
acoustic  waves  propagate  Independently  with  propagation  constants 
determined  by  Eqs.  (129)  and  (130),  respectively.  It  is  shown  later 
(Fig.  22)  that  the  acoustic  wave  is  attenuated  rapidly  in  this  region. 
Since  the  acoustic  wave  decays  to  a  very  small  fraction  of  its  ampli¬ 
tude  at  the  peak  within  a  very  short  distance  from  the  peak,  the  model 
is  shown  as  semi-infinite  in  the  figure.  The  temperature  profiles  are 
taken  to  be  linearly  increasing  also,  and  follow  either  Fig.  5(a)  or 
5(b)  in  form.  The  collision  frequency  is  assumed  to  be  a  constant. 

Boundary  Conditions 

First  consider  the  boundary  conditions  at  the  wall  (x  =  0) .  For 

convenience  (in  the  numerical  solution  and  the  approximate  solution) 

the  value  of  H  at  x  ■  0  is  considered  to  be  a  known  constant  (0.1  + 

2  Hz(0) 

jO.l  in  all  cases)  Hz(0).  Using  and  Poynting's  theorem  in 

the  dielectric,  and  assuming  that  an  area  of  1  in2  is  considered,  this 

value  of  H  (0)  corresponds  to  an  incident  electromagnetic  wave  with  a 
2 
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paver  of  from  6  w  to  24  nw  (for  R  =  1  or  B  *  0  above),  which  is  an 
easily  obtainable  power  level.  The  two  electromagnetic  boundary 
conditions  can  be  used  at  this  boundary  to  determine  the  amplitude  of 
the  incident  electromagnetic  wave  and  the  value  of  the  reflection 
coefficient.  Using  these  two  boundary  conditions  gives 


R 


n  .  cos  6,  H  (0)  -  E  (0) 
d  1  z  y 

n,  cos  0,  H  (0)  +  E  (0) 

n  4  2  V 


where  is  the  characteristic  impedance  of  the  dielectric,  and 


H 


v°> 

1  +  R 


(70) 


(71) 


where  Hq  is  the  amplitude  of  the  incident  electromagnetic  field.  The 
electric  field  at  the  boundary,  E^(0),  can  be  computed  from  the  and 
p  solutions  in  the  plasma  using  Eq.  (48).  The  acoustic  boundary  condi¬ 
tion  used  at  this  dielectric-plasma  interface  is  the  vanishing  of  the 
normal  component  of  the  electron  velocity  (V^).  Fran  Eq.  (49)  this 


results  in  a  specification  of  the  derivative  of  the  pressure  at  the 
boundary , 


iE 

dx 


eK  n  (0+) 
y  o 


x»Q 


we 


Hz(0+) 


(72) 
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where  nQ(0+)  is  the  electron  density  a.:  the  wall. 

Now  consider  the  boundary  conditions  at  the  distance  corresponding 
to  the  peak  in  electron  density.  The  two  electromagnetic  boundary 
conditions  are  the  same  for  either  of  the  two  plasma  profiles  consid¬ 
ered  in  Fig.  5.  The  acoustic  boundary  condition  will  be  different  for 
the  t  0  profiles.  For  the  profile  in  Fig.  5(a)  the  rigidity  boundary 
condition  is  used  which  imposes  a  zero  normal  component  of  the  electron 
velocity  at  t'lis  boundary.  This  implies  that  all  of  the  electrons 


4  i- 

caking  part  In  the  collective  motion  are  reflected  at  the  bourn  ary  and 
none  penetrate  Into  free  space.  Using  Eq.  (49),  the  derivative  of  the 
pressure  at  this  boixidary  Is  determined: 


where  n^(d  )  is  the  peak  electron  density  and  H^fd”)  is  the  amplitude 
of  the  magnetic  field  at  this  boundary.  For  the  profile  in  Fig.  5(b) 
the  dynamic  boundary  condition  la  used.  This  acoustic  boundary  con¬ 
dition  states  that  the  force  cm  one  side  of  the  bowdavy  must  balance 
the  force  on  the  other  side  of  the  boundary.  This  implies  that  the 
pressure  be  continuous  across  ths  boundary  except  at  an  interface 
between  a  compressible  medium  and  au  incompressible  medium.  The  kinetic 
boundary  condition  of  continuity  of  the  normal  component  of  the  elec¬ 
tron  velocity  ia  also  used. 


The  condition  of  continuity  of  at  the  outer  boundary  (x  -  d) 
determines  the  derivative  of  Uz  at  the  outer  botxidary.  For  the  pro¬ 
file  of  Fig.  5' -) ,  using  Eq.  (48)  tor  Ey  gives 

1  eKy  _  1 

~  Jwc(d-)  dx~  x>id_  +  uje(d-)A»  P^d  ^  juieQ  dx  xm£{+  '  ^ 


Now  assume  that  the  electromagnetic  field  for  x  >  d  has  only  an  outward 

-jK  x 
J  ox 

going  component  of  the  form  e  ,  where 


is  the  x  component  of  the  free  space  propagation  constant  for  electro¬ 
magnetic  waves.  Solving  Eq.  (74)  for  (dHz/dx)x<<d_,  gives 


*r  xcd.  -  j  s r  p<d'>  -  • 

It  la  shown  in  Appendix  C  that  for  the  peak  electron  densities  of  interest, 
the  second  term  on  the  left-hand  side  of  Eq.  (74)  is  completely  negligible 
compared  to  the  other  terms  in  the  equation.  Therefore,  to  a  very  good 


approximation 


dH 

sr  « -  V‘nKo*Vd+>  • 
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A  slightly  different  expression  results  for  the  derivative  of  Hz  at  the 
outer  boundary  for  the  profile  of  Fig.  5(b),  If  is  taken  to  be 
continuous  at  x  *  d  and  Eq.  (A8)  is  used  for  E^,  the  expression  for 
the  derivative  of  H  at  x  ■  d  becomes 
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Since  p  and  H  are  continuous  and  e  (d  )  ■  e  (d  )  at  x  *-  d, 
z  P  P 


*T  ,  ■  -'h,  1  -  (:8) 

x«d 

where  is  the  x  component  of  the  complex  propagation  constant  of  the 

electromagnetic  wave  in  the  homogeneous  plasma  region. 

Solutions  for  all  field  quantities  in  all  three  regions  can  be 

obtained  from  six  boundary  conditions,  three  at  each  boundary.  In 

the  numerical  solutions,  four  boundary  conditions  are  used  to  solve 


the  two  second  order  coupled  wave  equations  for  the  fields  in  the 
plasma  layer  for  the  profile  of  Fig.  5(a).  The  conditions  are: 
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(79) 


(90) 
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dx 
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J  p  c*  j' 


_z_2 —  a  m  ) 

VC  I  ' 


(31) 


{82) 


The  boundary  at  x  »  d  for  the  plasaa  profile  of  Fig.  5(b)  is  not 
associated  with  boundary  conditions  related  to  «a  interface  dividing 
two  different  aedin,  but  to  an  interface  where  the  fore  of  the  piasaa 
profile  changes  (the  undid*  recoins  die  sane).  One  extra  usknevn  is 
introduced  in  this  case — the  pressure  field  in  the  hoBogensoue  region 
which  did  not  exist  in  the  case  of  Fig.  5(a),  since  the  region  for 
X’  d  -jss  free  apace.  Inrcead  of  the  single  acoustic  betxidary  cccdi- 
ti'jr.,  tv o  Independent  acoustic  boundary  conditions  are  used  at  this 
interface.  The  first  oca  is  the  continuity  of  the  pressure  serose  the 
interface,  or 

p(d“)  -  p(d*>  (83) 

The  second  one  is  the  continuity  of  the  nomal  conpcoent  of  the 
ordered  electron  velocity  across  the  interface,  or 


y<r>  -  vx<d  )  (84) 

which  gives  (using  Eq.  (49)  and  the  continuity  of  H  at  x  *  d) 
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vbete  y  is  the  x  component  of  th*  complex  propagation  cons  tent  of 
t5»*  acoustic  «tvc  in  the  homogeneous  plasm*  region.  The  ordinary 
cleciwcwtgnetic  boundary  conditions  (E  and  H  continuous)  are  used  at 

y  * 

thu  *  -  d  boundary  for  profile  (b)  also.  The  boundary  conditions  used 
in  the  numerical  solution  for  profile  (b)  are; 
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P(<T)  -  p(d+)  . 
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Solutions  e t  Che  tvo  celled  «a*e  equations  sU!  b*  cbttloed  by 
npUdsg  the  wo  ar&mdrj  dlffecwtlAi  tqusUoof  with  tbcit  finite 
di  ffererce  app mineticas  -  ?h*  ore  equations  aw 
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vhere  the  principal  terse  are  or  Che  left-£and  aide  of  the  equations 
and  the  coaling  teres  are  on  the  right-hand  aide.  The  coefficients 
are  defined  as  follows: 
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Tbe  c<tf(lc>ecu  cf  tia  coupii-?  ^^Uoot 

tuiU  difference  ap^r^jciticn*  will  involve  -~a£p.:~x  oyasctities,  Cofc- 

oetelicn*  ievclrltj  the  complex  qusctitti*  aro  h*sd»<s?  *%asi  1-  tuifcg  the. 

cocplaa  atithmcic  capabi  line*  of  tha  FWT?Jtf  IV  * angcage  »  the  CSC 

6600  computer  The  v*ln  advantage  of  using  tCjo  tts»eric*l  solution  i* 

tfcJi  the  coupling  tetss  ere  retained  i«  the  fiek 

Since  the  u^aatlcns  m  ordinary  differential  ^satioes,  in  one 

iadcpecJsut  variable,  2,  the  first  step  is  to  divide  the  x-axis  into 

*  neJfcer  of  equally  spszei  points  s*  sbotm  In  Fig.  6.  The  step  sire 

Is  indicated  by  h.  In  the  finite  di ft* rone*  numerical  solution,  the 

variable  i  corresponds  to  discrete  points  os  the  x-axis.  Solutions 

for  H  and  p  will  consist  of  values  of  H  and  p  conouted.  at  x  values 
s  z 

corresponding  to  i,  from  the  finite  difference  approximations  to  the 
differential  equations.  The  fields  computed  at  the  discrete  points 
will  be  designated  as  and  p^. 

Letting  f  represent  either  field  quantity,  the  central  differ¬ 
ence  approximation  used  for  the  second  derivative  is  (Ref.  28) 


d2f  a  fi+l  "  2fi  *  ri-l 
dx2  h2 


(99) 


where  the  error  associated  with  this  approximation  is  equal  to 


where  |$|  <_  1  and  >U»  is  the  largest  value  of  the  fourth  derivative  of 
the  furctioo  in  the  region  of  interest.  The  finite  difference  approx¬ 


imation  used  for  the  first  derivative  is  the  central  difference  formula 
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USef.  2$) 


iL  ~  fi+i  ~  ?t-i 

dr  ~  2h 


(100) 


ti»*.  error  Mtodated  viii  this  approximation  U  equal  to 


■> 


where  j$j  «  1  and  Mj  Is  the  largest  value  of  the  third  derivative  in 
the  region  of  interest.  When  these  approximations  are  used  for  the 
derivatives  in  Eqs.  (91)  and  (92),  the  two  equations  become  (in  finite 
difference  form): 

K  -  2H  .  +  H  gt  ,,  Ih  -  H 

t(l-fl)  ti  z(l-lj  +  A  \  z(i+l)  ~  z(i-l) 

h2  t  *h 
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Pj+l  "  2pi  +  pi-l 

h2 


+  D 


+  EPi  «  ra* 


(102) 


The  solutions  for  E2i  and  new  require  the  solution  of  2(N+i)  linear 
algebraic  equations  where  N+l  is  the  number  of  points.  The  solution 
of  2 (N+l)  equations  is  required  since  and  must  be  determined 
from  Eqs.  (101)  and  (102)  at  a  total  of  N+l  points.  To  improve  the 
accuracy  of  the  tnxnerical  solution,  the  step  size  h  can  be  reduced 
(the  errors  are  proportional  to  h2) ;  however,  this  would  increase  N, 
requiring  the  solution  of  a  larger  number  of  algebraic  equations. 
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Ifc*  replacement  cf  tbe  different!.*)  equation*  by  the  ticiu  differ¬ 
ence  approximation* ,  £q*.  (lvl)  mo  (102) ,  is  ttzdard.  The  numerical 
technique  really  depend*  upon  the  method  chosen  for  the  solution  of  the 
resulting  algebraic  -equation*.  The  method  wed  hate  is  to  solve  for 
Btl  ana  fro*  the  finite  difference  equations.  Thin  glees 
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An  Iterative  method  for  solving  the  equations  is  used.  A  first  approx¬ 
imation  is  ured  to  calculate  a  second  approximation  which  in  turn  is 
used  to  calculate  a  third,  and  so  on.  Since  the  sane  formula  (Eqs. 

(103)  and  (104))  is  used  to  calculate  each  approximation,  the  Iterative 
process  is  said  to  be  stationery.  The  formulas  used  describe  a  point- 
method  since  the  next  approximation  at  one  point  is  explicitly  expressed 
in  terms  of  known  values  at  other  points.  The  iterative  procedure  is 
convergent  when  the  differences  between  the  exact  solution  and  the 
iterative  approximations  tend  to  zero  as  the  mmber  of  iterations 
increases.  The  Iterative  process  described  below  is  called  a  Jacobi 
iterative  method  (Ref.  29).  Initial  guesses  are  made  for  and  p^, 
across  the  entire  plasma  layer  (x  «*  0  to  x  *  d).  New  values  for 
and  p^  are  computed  at  each  point  using  Eqs.  (103)  and  (104)  and  the 
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initial  gaesm  for  and  p4 .  and  ;  are  again  computed  at  tads 

point  using  the  8^  and  p^  values  found  is  the  first  iteration.  The 
pro  cedar*  is  repeated  tntil  the  solution  converges  such  that  an  d  p. 

satisfy  Eqs.  (101)  and  (102)  at  all  points  across  the  larer.  When  this 
occurs,  the  and  p^  -confuted  fros  Eqs.  (103)  and  (104)  will  be 
or  changed  from  one  iteration  to  the  next,  this  will  never  happen,  of 
course,  because  of  die  errors  in  a  finite  difference  solution;  however, 
convergence  will  l>e  assumed  when  the  change  in  and  p  between 
succeaaive  iterations  be  cones  snail.  The  convergence  of  a  typical 
solution  is  Illustrated  in  Fig.  7  as  a  function  of  the  nibbtr  of  iter¬ 
ations.  The  convergence  of  the  p^  solution  is  such  faster,  as  shown 
in  Fig.  8,  A  typical  value  for  h  Is  d/50. 

Since  die  plasma  is  bounded,  the  field  solutions  axe  constrained 
to  obey  certain  conditions  at  die  boundaries  (see  the  previous  section. 
Boundary  Conditions).  These  are  included  In  the  numerical  solutions  by 
requiring  that  and  obey  these  conditions.  One  such  condition 
is  that  Hzo(i«0)  be  held  constant.  All  other  boundary  conditions 
involve  the  derivative  of  the  field  quantities.  The  derivatives  at 
the  boundaries  are  finite  difference  approximations  obtained  from  a 
Taylor  series  expansion  of  the  function  (Ref.  28)  at  the  appropriate 
boundary.  For  the  derivative  at  the  wall  (x  *  0) 
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where  the  number  subscripts  are  the  values  of  i.  For  the  far  boundary 
(x  -  d) 


the  derivative  values  will  i>e  used  to  solve  for  the  field  cuiJJti ties 
at  tint  botnoarieo.  Leasing 


where  Z}  to  Zu  are  coaplex  niabers.  The  value  of  H  at  x  «  0  is  fixed. 

z 

The  values  of  the  derivatives  of  the  functions  at  the  hoiBdaries  are 
used  tc  calculate  the  function  values  at  the  boundaries  (except,  of 
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couim,  for  the  raise  of  ^  it  i  *  3).  Equations  i  10S)  and  (106)  are 
used  to  calculate  tite  function  valued  at  tbc  beiaadaries  as  follows: 

— i;  -  Pi  *  *¥} 


Using  Eq.  U%)  st  a  «  d  fives 

i“-  '  rx-i  *  %-i 


„  **>  -  ".0.-2)  *  “.<*-)> 

"ril - -  • 

The  botadary  values  are  c^dattd  in  this  naoner  at  the  end  of  eadi 
iteration  and  are  need  for  calculations  in  the  sjcceeding  Iteration. 

Coefficient*  of  the  Coupled  Equations 

The  coepltz  coefficients  of  tbc  coupled  wave  equations  (61)  and 
(62)  determine  the  forms  of  the  solutions  for  or  p.  The  importance 
of  die  magnitude*  of  the  coefficients  is  illustrated  in  the  neat  section. 
By  woririnf  with  the  coefficients  for  various  plasma  profiles  and 
operating  frequencies ,  one  finds  that  certain  coefficients  are  much 
more  Important  than  others.  The  coupling  coefficient  F  is  perhaps  the 
most  significant  coefficient.  When  the  easet  form  of  the  continuity 
equation  (Eq.  (10))  is  used,  coefficient  F  is  given  by 


F 


eK 


—I 


we 
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(107) 


This  coupling  coefficient  as  given  by  Bureau  (Kef.  16)  includes  only 


tnt  tiist  i«ts.  (fee-  baclttt  vJttuK  <»f  the  ccw-fficfent  F  given  by  !«. 

(>''•?  ana  by  turMB  <*ef.  Ife)  *ne  plotted  i«  Fig.  9.  tz  is 

ats«-7i74  that  Utr  r  iDiludisi  tue  adrfitLc&al  lers  is  larger,  especially 

at  t  t  Mtiitr  clectmi  deulries.  >ince,  In  the  ishoH|(ato»  region . 

the  f.n*iS  JK  is  essentially  prcfortttecl  to  f,  the  ugnitadt  of  F  could 

m*c  a  large  difference  la  the  pressure  spliisde.  To*  other  coupling 

coefficient  i*  plotted  in  Fig.  10.  loth  coupling  coefficients  reach 

a  tt_LsIauE  vsiue  *t  the  sane  dlstaace  into  the  plan**.  The  peek  ie  the 

curve  representing  coefficient  F  is  caused  by  *  stolsai  in  c  St  3 

P 

point  where  .  »  , ;  i.e.,  *n  x  rslue  where  the  electron  density  a 

p  -  o 

results  ie  a  value  of  b 

P 


equal  to  the  operating  frequency,  Since  «  Is  assured  to  be  linearly 
increasing,  die  x  dependence  of  F  is  given  by  oo/£p;  therefore,  the 
coupling  coefficient  is  directly  proportional  to  the  density  and 
inversely  proportioned,  to  the  relative  dielectric  constant  of  the 
piasRa  given  by 


c 

P 
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for  a  lossy  plasma  and 
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Figure  10.  Coupling  Coefficient  C 
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far  a  Iocs  Less  pljssi.  For  a  lossless  plasma  c  *0  for  u  »  .  and  F 

P  P 

vould  be  infinite;  however,  c  hts  s  Minis  us  for  -  *  .  is  *  lossv 

P  P 

plns»  and  F  has  •  peek  value  at  the  distance  z  where  *  v.  The 
peak  will  be  less  sharp  for  larger  collision  frequencies. 

The  other  coefficients  are  given  in  Appendix  B  for  a  typical 
plana  profile  aid  operating  frequency. 

Approximate  Soli tloc 


In  determining  an  approximate  solution,  the  relative  magnitudes 
of  the  terms  in  the  equations  are  considered.  If  some  teems  In  the 
equations  are  much  smaller  than  others,  these  can  be  neglected  and  the 
equations  simplified.  The  procedure  la  described  belm: 

(1)  Replace  a  second  derivative  term  by  the  dependent  variable 
in  die  derivative  divided  by  the  square  of  the  characteristic  length 
(the  plazas  thick*****,  d) . 

(2)  Replace  a  first  derivative  tern  by  die  dependent  variable 
la  the  derivative  divided  by  the  characteristic  length. 

(3)  Determine  the  relative  magnitude  of  the  resulting  teres 
from  the  coefficient  magnitudes. 

For  the  two  coupled  differential  equations  the  first  two  steps 
result  in: 


(108) 


(109) 
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Fur  an  electron  density  profile  which  increases  linearly  froe  10! ' 
elecirons/cs  to  JO^  clectrons/ca*  In  2  as,  it  is  shown  in  Appendix 
A  chat  for  certain  frequencies  the  above  equations  becone 


and 


(1J0) 


Ep  -  rdz  Uii) 

with  the  other  tents  being  such  smaller.  Substituting  for  p  in  Eq. 
(110)  fro*  Eq .  (Ill)  gives 


or 


(112) 


(113) 


It  is  also  shown  in  Appendix  A  that  CF ft  is  negligible  compared  to 
the  other  tens  in  the  parenthesis  in  Eq.  (113).  Equations  (108)  and 


(109)  thus  becone 


and 


Reintroducing  the  derivative  citation,  these  equations  bectsae 


d2H 

z 

dx2 


0 


U24) 
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•(  l k  • 


(115) 


This  procedure  appears  similar  to  a  successive  approximation  technique; 
however,  it  should  be  coted  that  the  coupling  terms  were  neglected 
ooly  when  their  magnitude  was  much  smaller  than  the  other  terms  in  the 
equation.  Is  the  equation  the  coupling  term  was  negligible;  however, 
in  the  p  equation,  the  coupling  term  was  a  dominant  factor. 

Equation  (114)  can  be  solved  for  H£.  The  pressure  p  can  then  be 
fovnd  from  £q.  (115).  First  define 


dx  * 


where 


♦  AG  *  0  , 


i  dc 

P 


dG  _  / 1_  djA 

e  Vp 


dx  ■  0  . 


Integrating  both  sides  of  the  above  equation  gives 


in  G  ■  in  e  +  in  C  , 

P 


G  -  C;  ep  , 
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drt 

— -  «  C  * 
ax  P 

where  c;  is  the  first  constant  of  integration. 

Integrating  both  sides  of  the  above  equation  gives 

s 

H  •  C  -.  I  [  dx  +  C  -. 

*  ‘J  P 

o 

where  £,  is  tire  second  constant  of  integration. 

Since  H  (x<*0}  «  H  (0) 

2  2 

X 

Ht  -  Cl  J  Cp  dx  +  lyo)  .  (life) 

o 


Using  che  rigidity  boundary  condition  at  x  •  d.  It  is  shown  by  Eq, 
(76)  that 


!s»d 


-  -It  (dT)  & 

J  p  ax  z 


H  (<^> 


Using  this  condition  gives 


Ci 


(117) 


where  e^(x)  is  the  relative  permittivity  (complex)  of  the  plasma  evalu¬ 
ated  at  x,  and  Kqx  is  given  by  Eq.  (75). 

The  approximate  solution  is  not  always  accurate  for  the  plasmas 
of  interest;  however,  the  approximate  p  solution  can  be  used  for  all 
profiles,  since  coefficient  E  is  very  large  (see  Appendix  B).  The 
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approximate  solutions  «rs  compared  with  ntxaericai  solutions  in  figs. 

11  and  12.  Io  cost  casts,  the  initial  guesses  for  a  nuaerical  solu¬ 
tion  were  the  approximate  solutions  described  above. 

Solutions  of  the  Coupled  Wave  Equations 

The  Meaning  of  the  and  p  solutions  will  be  discussed  before 
the  results  are  presented.  The  field  solutions  are  valid  solutions 
of  the  coupled  wave  equations;  however,  the  curve*  describing  the 
pressure  field  as  a  function  of  a  shew  only  a  part  of  the  solution. 

If  the  coupling  ten  on  the  right-hand  side  of  Eq.  (61)  is  neglected, 
the  pressure  field  solution  has  two  parts.  One  part  consists  of  the 
mathematical  representation  of  two  plane  waves  traveling  in  opposite 
directions  (attenuated  by  the  same  anoint  as  a  function  of  x)  and 
the  other  part  consists  of  a  pressure  amplitude  with  a  different  x- 
depeodenos.  It  will  be  shown  that  the  results  presented  In  this 
chapter  give  only  the  second  part  ef  the  pressure  solution. 

As  discussed  in  the  section  on  the  approximate  solution,  the 
coupling  (or  influence)  of  the  tern  on  the  right-hand  side  of  Eq.  (61) 
Is  negligible  for  nany  of  the  profiles  and  frequencies  of  Interest. 

The  complete  analytical  solution  for  p  will  be  obtained  by  making  some 
rather  broad  assumptions,  hence  the  error  in  neglecting  the  coupling 
term  in  Eq.  (61)  will  not  present  serious  problems.  The  equations  to 
be  considered  are: 

d2H  dH 

- *  +  A  -tt-  +  BH  =  0  (118) 

dx2  dX  2 


and 
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Figure  11.  Comparison  of  Approximate  Analytical  Solution 
anci  Numerical  Solution,  Pressure  Field 


Figure  12.  Comparison  of  Approximate  Analytical  Solution 
and  Nunerical  Solution,  Magnetic  Field 
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+  «  (119) 

d*  * 

wJ.->re  g  *  FF^.  The  sagcetic  field  H  can  be  determined  from  the  firat 
equalise.  The  second  equation  la  now  considered  in  detail.  This  is 
an  inhomogeneous ,  second  order,  ordinary  differential  aquation.  The 
solution  will  consist  of  two  parts:  the  homogeneous  solution  and  the 
particular  solution.  If  D  and  E  are  constants  (whidh  is  not  a  had 
asaxpiioQ  for  a  linearly  increasing  electron  density  profile),  the 
homogeneous  solution  is 

y  *  yt* 

ph  -  ee  1  ♦  be  ‘  ,  (120) 


where  y  and  *  are  the  two  roots  of  the  equation  a2  ♦  Os  ♦  E  •  0  and 

J  2 

.  mt  b  ...  «—««..  na. 

is  due  to  die  large  collision  frequency  in  the  plai 
The  particular  solution  is  calculated  aa  follows: 


(see  Eq.  (97). 
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(121) 


where  W(6)  can  be  determined  by  calculating  die  Vronski  an 


W(x) 
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(Tr%): 


and  the  Vronski an  is  independent  of  x.  Tbe  particular  solution  is 
disc  fond  to  be 


where  it  on  assumed  that  she  product  FH^  »  g  was  constant,  allowing 
the  integrate of  Eq.  (121)  to  be  perforned  easily.  The  complete 
solution  is  jnt  the  sis  of  the  homogeneous  solution  end  die  partic¬ 
ular  solution 


which  can  be  rewritten 

>1*  Y2x 

P  -  A'e  +  B'e  +  */E  ,  (122) 

where  YjY2  ■  E»  and  A'  *  a  -  g/2E  and  B*  ■  b  -  g/2E  are  defined  as  the 
new  constants  of  rhe  homogeneous  solution  (since  g  and  E  were  considered 
to  be  constant).  The  constants  A1  and  B'  can  be  deterained  from  the 
two  boundary  conditions  and  are  found  to  be  functions  of  both  and 
dHg/dx  at  both  bomdaries.  The  particular  solution  has  no  undetermined 


70 


coefficient*  bet  depends  cm  the  "forcing  function**  g.  The  particular 
solution  also  is  a  fiaetioe  of  H  since  «  *  FH  . 

The  niserical  solution  yields  pressure  solutions  with  approxi¬ 
mate  ly  the  foilobiLg  dependence 

P  -(  |  )Ht  -  g/E  ,  (123) 

and  diis  is  just  the  first  pare  of  the  particular  solution.  In  Che 
solution,  the  other  terns  Kust  be  negligible  since  the  a  dependence 
is  given  by  g/E  alone  and  the  coup  let e  solution  is  given  by  the  cub  of 
the  terns.  It  is  concluded  that  the  "forced  solution'"  (s  tern  in  the 
particular  solution)  is  such  larger  than  the  terns  in  the  solution 
characterizing  the  acoustic  wave  propagation  in  the  plssna. 

It  was  noted  that  changing  the  acoustic  boimdary  conditions 
(e.g.t  setting  p  •  0  at  *  -  d)  had  very  little  effect  on  the  solution 
at  interior  points.  This  should  be  expected  since  the  particular 
solution  is  not  affected  by  changing  the  const sots  A'  and  B'  which 
depend  on  the  boundary  conditions. 

Even  though  the  two  plssna  models  considered  are  sinple,  there 
are  still  eight  different  paraseters  which  can  be  changed  to  affect 
the  solutions.  These  Include:  the  electron  density  at  the  wall, 

Che  peak  electron  density,  the  tenperature  at  the  wall,  the  tenpera- 
ture  at  the  distance  corresponding  to  the  peak  electron  density,  the 
collision  frequency,  the  distance  to  the  peak,  the  operating  frequency 
and  the  angle  of  incidence.  For  the  results  presented  in  Figs.  13  to 
22,  the  angle  of  incidence  was  fixed  at  30°. 


?i 

In  Fig*.  13  to  18  the  profile  : f  Fig.  5(a)  is  wei  Th«.  vail  «a<J 
peak  electron  acuities,  the  wall  and  peak  temperatures,  sad  the 
collision  frequency  are  indicated  on  each  figure.  Solutions  corres¬ 
ponding  to  different  operating  frequencies  are  labelled  with  the 
particular  frequency.  The  a ague tic  field  solution  to?  the  first  desaity 
profile  is  given  in  Fig.  14  for  three  differeut  operating  frequencies. 
Since  the  scoustl-  aodes  generated  in  the  plasma  are  of  priaary  concern, 
this  is  the  only  aagnetic  fle*d  solution  presented.  Figure  15  gives 
the  pressure  field  for  the  sane  wall  and  peak  electron  densities  as 
Fig.  13;  however,  the  distance  to  the  peck  in  Fig.  15  la  only  0.5  hi 
as  coapared  to  1.0  aa  in  Fig.  13.  This  iaplies  that  the  gradient  of 
the  electron  density  in  Fig.  IS  Is  larger  than  that  in  Fig.  13.  Fig¬ 
ure  16  gives  pressure  solutions  for  a  low  electron  density  at  the  wall 
(1010  cu"3) ,  while  Figs.  17  and  18  give  pressure  solutions  for  a  high 
peak  electron  density  (10* '*  ca~5) .  As  given  by  Eq.  (123),  the  x- 
dependence  of  p  is  just  FH^/E.  Since  E  is  almost  constant  across  the 
plssas  layer  (see  Appendix  B)  and  H  is  slowly  varying  with  no  pro- 
nouaced  peaks,  the  peaks  in  the  curves  representing  the  pressure  solu¬ 
tions  are  due  to  F.  This  peak  in  F  can  be  seen  in  Fig.  9.  The  peak  oc¬ 
curs  at  the  distance  where  the  electron  plasaa  frequency  (determined  by 
die  electron  density)  is  equal  to  the  operating  frequency.  The  magni¬ 
tudes  of  the  pressure  solutions  for  the  different  density  profiles  are 
not  significantly  different;  however,  comparing  the  30  GHz  curve  in  Fig. 
13  to  the  30  GHz  curve  in  Fig.  15,  it  is  noted  that  the  pressure  magni¬ 
tude  for  the  larger  electron  density  case  (Fig.  15)  is  slightly  larger. 
The  significant  result  of  the  pressure  solutions  given  in  Figs.  13  to 
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Figure  16.  Pressure  Field  for  1010  -  1013  cm”3 
Density  Profile 
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Figure  17.  Pressure  Field  for  1011 
Density  Profile 
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18,  as  stated  earlier,  is  that  the  curve  clearly  describes  a  variation 
given  bv  Eq.  (123)  with  no  x- variation  of  the  type  given  by  the  two 
exponential  terms  in  Eq.  (122)  observable.  The  influence  of  these 
terms  (describing  acoustic  wave  propagation)  on  the  "forced  solution," 
Eq.  (123),  must  be  negligible.  From  Figs.  13  to  18,  it  is  observed 
that  the  amplitudes  of  these  terms  must  be  of  the  order  of  10“8  newtons/ 
meter  or  smaller.  The  amplitudes  of  these  terms  are  discussed 
further  in  Chapter  V. 

In  Fig.  19,  the  pressure  field  solutions  are  given  for  various 

temperature  profiles  for  a  1011  -  10 13  cm“3  electron  density  profile. 

The  pressure  increases  as  the  plasma  temperature  is  Increased. 

In  Fig.  20  the  profile  of  Fig.  5(b)  is  used.  Comparing  the  30 

GHz  curve  in  Figs.  13  and  20,  it  is  noted  that  the  pressure  amplitude 

for  profile  (b)  is  smaller  than  the  amplitude  for  profile  (a).  This  is 

due  to  the  fact  that  Hz  is  smaller  for  (b) .  The  same  x-dependence 

(FH  /E)  is  present  for  the  pressure  solutions  in  both  profiles.  The 
z 

same  explanation  holds  for  Fig.  21. 

The  attenuation  of  a  pressure  wave  propagating  in  the  positive 
x  direction  from  x  ■  d  in  the  homogeneous  region  of  profile  (b)  was 
calculated  by  determining  the  real  part  of  the  x  component  of  the 
propagation  constant  given  by  Eq.  (130).  The  pressure  amplitude  was 
normalized  at  x  ■  d  and  the  attenuation  is  described  by 

I  -  e'“<’l'd)  . 

where  a  is  the  attenuation  factor  from  Eq.  (130).  The  wave  is  attenu¬ 
ated  by  a  factor  of  lO-4  within  a  distance  of  only  0.02  mm  from  the 
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peak  (see  Fig.  22).  The  calculations  were  made  at  low  frequencies 
(1  GHz  and  10  GHz)  where  the  attenuation  is  less  severe.  Higher 
operating  frequencies  would  result  in  an  even  more  rapid  damping  of 
the  wave . 


CHAPTER  V 


CONVERSION  EFFICIENCIES  AT  A  DIELECTRIC-PLASMA  INTERFACE 
Conversion  of  an  Electromagnetic  Wave  into  an  Acoustic  Wave 

Ic  has  been  shown  chat  the  amplitudes  of  waves  associated  with  the 
propagating  electron  acoustic  mode  are  of  the  order  of  i0~9  newtons/ 
meter^  or  smaller  even  in  a  plasma  region  where  there  exists  coupling 
between  the  electromagnetic  wave  and  the  acoustic  wave.  An  estimate 
of  the  wave  amplitudes  is  made  in  the  section  Power  Conversion  at  the 
end  of  this  chapter.  Even  though  it  appears  that  the  amplitude  of  the 
reflected  acoustic  wave  is  so  small  as  to  be  undetectable,  the  conver¬ 
sion  efficiencies  will  be  investigated,  since  very  few  numerical 
results  seem  to  be  available  even  for  the  semi-infinite  plasma  case 
of  this  chapter.  For  the  semi-infinite  plasma,  the  pressure  solution 
simply  consists  of  one  of  the  exponential  terms  in  the  homogeneous 
solution  of  Chapter  IV  which  corresponds  to  propagation  in  the  posi¬ 
tive  x  direction. 

An  electron  acoustic  wave  can  be  generated  at  the  boundary  between 
a  dielectric  material  and  a  plasma  by  an  electromagnetic  wave  which  is 
vertically  polarized  and  obliquely  incident  from  the  dielectric  mate¬ 
rial  (Ref.  21).  It  was  shown  in  Chapter  III  that  a  vertically  polar¬ 
ized  electromagnetic  wave  incident  on  an  inhomogeneous  plasma  generates 
a  disturbance  described  by  two  differential  equations  implying  a 
coupling  between  an  electromagnetic  mode  and  a  plasma  mode  because  of 
the  inhomogeneity  of  the  mediun.  The  modes  are  also  coupled  at  a 
boundary,  but  the  coupling  in  this  case  is  through  the  boundary 
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conditions  which  are  both  electromagnetic  and  acoustic.  The  boundary 
conditions  used  are  the  Continuity  of  the  tangential  components  of  the 
electric  and  magnetic  fields  (Ref.  30)  and  the  vanishing  of  the  normal 
component  of  the  electron  velocity  (Refs.  12,  13,  and  22). 


H  (O')  -  H  (0+) 

z  z 

(124) 

E  (0”)  -  E  (0+) 

(125) 

y  y 

V  -  0 

X 

(126) 

The  acoustic  boundary  condition,  Eq.  (126),  implies  that  none  of  the 
electrons  penetrates  the  boundary  and  is  sometimes  called  the  "rigidity 
boundary  condition."  This  acoustic  boundary  condition  is  considered 
to  be  valid  at  a  dielectric  plasma  interface  (Refs.  12,  13,  21,  and 
31);  however,  at  the  boundary  between  a  metal  and  a  plasma,  a  different 
acoustic  boundary  condition  must  be  used  (Ref.  31)  which  specifies  the 
ratio  of  the  normal  component  of  the  electron  velocity  to  the  pressure 
and  is  referred  to  as  the  "impedance  boundary  condition." 

In  this  chapter,  only  homogeneous  plasmas  are  considered  so  that 
the  spatial  dependence  of  the  fields  can  be  explicitly  written.  Plane 
wave  propagation  is  considered  with  harmonic  time  variation  assumed. 

The  objective  is  to  determine  nunerically  the  effects  of  plasma  param¬ 
eters  (electron  density,  temperature,  and  collision  frequency)  and  the 
incident  electromagnetic  wave  (frequency  and  angle  of  incidence)  upon 
the  conversion  efficiencies  at  the  boundary. 

The  conversion  of  an  electromagnetic  wave  into  a  pressure  wave  at 
a  dielectric-plasma  half-space  interface  is  shown  in  Fig.  23.  The 
electromagnetic  wave  is  vertically  polarized  and  incident  at  some 
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(132) 
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where  the  y  dependence  is  the  same  as  that  for  Eqs.  (52)  and  (53). 

The  value  of  the  pressure  wave  at  the  boundary  Is  p(0)  and  H  Is  the 

zo 

value  of  the  transmitted  magnetic  field  at  the  boundary. 

The  electromagnetic  field  In  the  dielectric  is  taken  to  be  a 
superposition  of  an  incident  plane  wave  and  a  reflected  plane  wave 


and  is  expressed  in  terms  of  the  magnetic  field  as 


/ 

h  «=  h  (< 

z  o  \ 


"YdxX  ^  „  Ydx: 
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where  is  the  amplitude  of  the  incident  wave  and  R  is  defined  as  the 
reflection  coefficient.  The  x  component  of  the  propagation  constant 
in  the  dielectric  is  defined  as  y »x  and  is  given  by 


Yj  »  K2  -  (o2u  e  K  , 
'  dx  y  o  o  d 


(134) 


y  *  luiYy  c  K,  cos  0,  . 
1  dx  J  o  o  d  1 


(135) 


The  ratio  p(0)/H  will  be  determined  using  the  three  boundary 
o 

conditions,  Equ.  (124),  (125),  and  (126). 

From  Eq.  (124)  (Hz  continuous  at  the  boundary  x  =  0) , 


H  +HR-H 
0  0  zo 


(136) 


From  Eq.  (125)  (E^  continuous  at  the  boundary  x  =  0) ,  and  using 
Ey  =  -l/jwe^  ^dHz/dxj  in  the  dielectric  along  with  Eq.  (133)  and  using 
Eq.  (48)  for  E^  in  the  plasma  along  with  Eqs.  (131)  and  (132), 
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From  Eq.  (126)  (V^  *  0  at  the  boundary  x  *  0) ,  and  using  Eq.  (49) 

for  V 
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Equations  (136),  (137),  and  (133)  can  be  used  to  solve  for  the 
ratio  p(0)/HQ,  which  gives  the  magnitude  of  the  generated  pressure 
wave  at  the  bowdary  relative  to  the  amplitude  oi  the  Incident  electro 
magnetic  wave.  The  result  is 
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Equation  (139)  will  be  used  to  investigate  tne  conversion  of  a  trans¬ 
verse  electromagnetic  wave  into  a  pressure  wave  at  a  plane  dielectric- 
plasma  interface.  A  computer  program  was  written  to  solve  Eq.  (139) 
for  various  plasma  parameters  and  incident  electromagnetic  wave 
properties.  The  absolute  value  of  Eq.  (139)  is  used.  For  simplicity, 
this  value  will  be  designated  as  | p /H |  in  the  following  graphs. 

The  electromagnetic-to-acoustic  conversion  efficiency  is  plotted 
as  a  function  of  electron  density  in  Fig.  24.  The  temperature  and 
collision  frequency  will  vary  for  different  electron  densities.  Col¬ 
lision  frequencies  and  temperatures  for  different  electron  densities 
are  shown  in  Table  I.  These  values  were  taken  from  case  studies  made 
of  the  plasma  sheath  surrounding  a  hypersonic  reentry  vehicle  u3ing  a 
computer  code  developed  by  McDonnell  Douglas  Corporation  (Ref.  32). 


ELECTRON  OENSITY 
(ElECTRONS/m*) 


Figure  24.  Electromagnetic-to-Acoustic  Conversion 
Efficiency  versus  Electron  Density 


Table  I 


COLLISION  FREQUENCIES  AND  TEMPERATURES 


FOR  GIVEN 

ELECTRON  DENSITIES 

(Ref.  32) 

Electron 

Collision 

density 

frequency 

Temperature 

(cm"  ) 

(sec"  ) 

(°K) 

1.6  x  1011 

10n 

2000 

3.2  x  1011 

10n 

2000 

6.4  x  1011 

10n 

2000 

1.3  x  1012 

10 11 

3000 

2.6  x  1012 

101! 

3000 

5.1  x  1012 

101* 

3000 

1.0  x  ID13 

1011 

4000 

2.0  x  1013 

2  x  10’* 1 

4000 

4.1  x  1013 

2  x  10 31 

4000 

6.2  x  1013 

4  x  1011 

4000 

The  program  treats  the  invlscld  flow  field,  ablation  mass  less,  and 
turbulent  boundary  layer  surrounding  a  spherically  blunted  conical 
body.  The  electron  density,  temperature,  and  collision  frequency  in 
the  turbulent  boundary  layer  are  of  primary  concern  in  this  report. 
The  case  studies  were  made  by  the  Fuzing  Environment  Brandi  (WLEE)  of 
the  Air  Force  Weapons  Laboratory  using  the  computer  code  for  several 
advanced  reentry  vehicles  with  two  different  types  of  heatshield 
materials  and  various  reentry  trajectories.  The  significant  result 
of  Fig.  24  is  that  the  electromagnetic-to-acoustic  conversion  effi¬ 
ciency  increases  monotonically,  with  increasing  electron  density. 
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fhis  tallies  chat  a  sort  dense  flwu  would  result  la  a  better  coupling 
of  elect  wag— tic  energy  to  acoustic  energy  at  the  boucdary.  The 
dashed  portion  of  the  40  (31s  «r»  is  Fig.  24  Indicates  the  region 
where  the  acoustic  vases  axe  Lsidsi  dsgei  sines  lp2  in  this  region. 

Figure  25  flume  the  elect  tomagmetlc-to  acoustic  const r* loo  effi¬ 
ciency  as  a  function  of  tbs  plassis  temperature.  It  is  noted  that  the 
conversion  efficiency  increases  as  the  plasms  temperature  in  increased. 
This  is  to  be  expected  sines  it  is  the  plasms  tcnpsrstuxe  upon  which 
the  acoustic  vases  depaad  as  s  propagating  mode  (see  Eq.  (33)).  The 
numerical  solutions  show  that  the  "forced  solution"  of  the  pressure 
field  also  Increases  with  increasing  temperature  as  illustrated  in 
Fig.  19. 

Figures  26,  27,  and  28  iilustrata  the  fact  that  the  electromagnetic 
to  acoustic  conversion  efficiency  has  a  maximal  value  at  a  particular 
angle  of  incidence.  For  on  electron  density  of  10 1 3  cn-3  and  a 
frequency  of  40  GHz  the  efficiency  ban  a  peak  at  approximately  25" 

(Fig.  26).  For  a  frequency  of  31.25  GHZ  the  piak  still  occurs  at  25* 
(Fig.  27);  however,  Che  efficiency  aagnltude  is  larger  Ly  an  order  of 
magnitude  due  to  tha  frequency  dependence  of  die  efficiency  (see  Figs. 

29  and  30).  When  the  electron  density  is  Increased  to  l.'*4  an-3,  the 
peak  occurs  at  45*  (Fig.  28). 

Figures  29  and  30  illustrate  the  frequency  dependence  of  the 
electromagnetic-to-acoustic  couversion  efficiency.  The  efficiency  is 
plotted  as  a  function  of  the  dimensionless  parameter  w/w^,  where 
is  the  electron  plasma  frequency  detemined  by  the  electron  density 


40  GHz 
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Figure  29.  E le ct rocugne ti c-to-Acous tic  Conversion  Efficiency 
versus  <u/u _),  (n  «  10 12  cn"a) 
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from  Eq.  (30).  The  curves  exhibit  a  minimus  at  a  frequency  slightly 
above  the  plana  frequency  and  a  steady  decrease  at  frequencies  greater 
than  about  twice  the  plana  frequency.  The  physical  explanation  for 
this  behavior  lies  in  the  interpretation  of  the  electron  plans  fre¬ 
quency.  At  frequencies  below  the  plasaa  frequency,  the  electrons  can 
respond  to  the  tine  variation  of  the  electrooagnetic  wave  and  at  fre¬ 
quencies  above  the  plasna  frequencies  the  electrons  cannot  follow  the 
field.  Since  the  pressure  wave  describes  a  collective  notion  of  the 
electrons,  a  larger  response  is  expected  at  frequencies  below  the 
plasna  frequencies.  Fron  Eq.  (139)  it  is  seen  that  the  efficiency  is 
directly  proportional  to  the  relative  dielectric  constant  of  the  plasna, 

c  .  For  e  collisionleaa  plasna,  e  is  zero  for  w  «  u  and  the  effi- 
P  P  P 

ciency  would  be  zero.  For  a  plasna  with  a  large  collision  frequency 

such  as  the  one  considered  here,  c  has  a  minimus  for  u  ■  w  and  does 

P  P 

not  vanish.  The  dashed  portion  of  the  curve  in  Fig.  29  indicates  that 
Landau  damping  would  destroy  any  collective  motion  of  the  electrons  in 
this  frequency  range.  The  efficiency  is  an  order  of  magnitude  smaller 
for  an  electron  density  of  1012  cm“3  (Fig.  29)  as  compared  to  an  elec¬ 
tron  density  of  1014  cnT3  (Fig.  30). 

The  electromagnetic-to-acoustlc  conversion  efficiency  is  plotted 
as  a  function  of  the  dimensionless  parameter  v/u>  in  Fig.  31.  The 
effect  of  a  large  collision  frequency  is  to  reduce  the  amplitude  of 
any  organized  motion  of  the  electrons  considered  as  a  continuun.  The 
effect  on  the  conversion  efficiency  is  therefore  a  reduction  in  magni¬ 
tude  of  the  conversion  efficiency.  The  reduction  is  fairly  small  and 
constant  for  collision  frequencies  smaller  than  the  operating  (radian) 


Electrooagnetlc-to-Acouatic  Conversion  Efficiency  versus  (v/w) 
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frequency  (Fig.  31);  howaver,  Che  reduction  in  the  efficiency  magni¬ 
tude  is  greater  for  collision  frequencies  larger  than  the  operating 
(radian)  frequency  as  illustrated  in  Fig.  31. 


Conversion  of  an  Acoustic  Wave  into  an  Electnanagnetic  Wave 

An  electron  acoustic  wave  incident  upon  a  boundary  between  a 
plasma  and  a  dielectric  can  excite  an  electromagnetic  wave  in  the 
dielectric.  This  is  a  consequence  of  the  boundary  conditions  being 
both  electromagnetic  and  acoustic  in  nature,  Eqs.  (124)  to  (126).  The 
conversion  process  is  illustrated  in  Fig.  32.  In  general,  there  will 
be  both  transmitted  and  reflected  electromagnetic  waves  excited  at  the 
boundary.  As  shown  in  Chapter  III,  a  pressure  wave  such  as  the  one 
shown  in  Fig.  24  can  be  excited  only  by  a  vertically  polarized  electro* 
magnetic  wave,  therefore,  a  pressure  wave  of  this  type  will  excite 
only  a  vertically  polarized  electromagnetic  wave. 

The  plasma  is  considered  homogeneous  so  that  expressions  for 
the  fields  can  be  written  explicitly.  There  is  one  other  aasuaption 
which  must  be  made  in  order  to  obtain  nimerical  results.  The  y  depend¬ 
ence  of  all  field  quantities  (Including  the  incident  pressure  wave)  is 
as  sided  to  be 


(140) 


where  K 

y 


is  given  by 


■  ui/y  e  K,  sin  0i  . 
o  o  a 


K 

y 


(141) 
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The  value  of  will  determine  the  angle  9  at  which  the  pressure  wave 

la  incident  upon  the  b oiai da ry ,  since  K  can  be  determined  from  y^x 

(K  is  the  isuglaary  part  of  y  )  and 
px  px 

%  ■  ■“"[VpJ  • 

The  spatial  variations  of  the  field  quantities  are  determined  in 
the  same  manner  as  in  the  last  section.  The  electromagnetic  wave 
reflected  in  the  positive  x  direction  is  described  by  the  magnetic 
field  strength  in  the  plasma 
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where  Hzr  is  the  amplitude  of  the  reflected  magnetic  field  strength 
at  x  -  0  and  yhx  is  defined  by  Eq.  (129).  H  has  only  a  z  component, 
for  a  verticrJLly  polarized  electromagnetic  plane  wave.  The  pressure 
wave  in  the  plasma  is  given  by 
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where  p^  is  the  amplitude  of  the  incident  pressure  wave  at  x  *•  0  and 

p  i"  the  amplitude  of  the  reflected  pressure  wave  at  x  *  0  an!  y 
r  px 

is  defined  by  Eq.  (130).  The  transmitted  electromagnetic  wave  in  the 
dielectric  is  described  by 
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where  H  is  the  amplitude  of  the  transmitted  magnetic  field  strength 
at  x  «  0  *id  Is  defined  by  Eq.  (135). 

Applying  the  botadary  condition  froa  Eq.  (124)  gives 


U  «  H 
zt  zr 


(142) 


Applwins  the  boundary  condition  from  Eq.  (125)  gives 


eKy  ’ 

uc Ait  (pi 


(143) 


litis  completes  the  two  electromagnetic  boundary  conditions. 

Applying  the  acoustic  boundary  condition  froa  Eq.  (126)  gives 


eK  e 

dk  "zr  -  TiST  (v  pi  '  V  pr)  ■  0  ’  (li 

where  Eq.  (49)  has  been  used  for  V^. 

Equations  (142),  (143) t  and  (144)  can  be  t  .ved  (or  the  ratio 
H^/pj  which  gives  the  ratio  of  the  aaplitude  of  the  electromagnetic 
field  excited  in  the  dielectric  to  the  aaplitude  of  the  incident 
pressure  wave.  The  result  is 


(145) 


Equation  (145)  will  be  'jped  to  calculate  the  conversion  of  a 
pressure  wave  into  an  electromagnetic  wave  at  a  plane  dielectric- 
plasma  interface  for  varying  plasma  parameters  A  computer  program 


10S 

warn  written  to  solve  Eq.  (145).  The  absolute  value  of  the  ratio  in 
Eq.  (145)  is  used.  This  value  will  be  designated  as  iH/pj  in  the 
following  figures. 

The  acoustic-to-eiectroaagnetic  conversion  efiicicncy  is  plotted 
as  a  fraction  of  electron  density  for  an  operating  frequency  of  2  GHz 
in  Fig.  33  and  for  40  GHz  in  Fig.  34.  The  collision  frequencies  and 
temperature*  corresponding  to  the  electron  densities  at  which  the 
efficiency  was  computed  are  given  in  Table  1.  The  conversion  effi¬ 
ciency  is  larger  for  lower  electron  densities.  This  is  desirable, 
since  the  conversion  process  (from  a  pressure  wave  into  an  electro¬ 
magnetic  wave)  will  occur  at  the  dielectric-plasma  interface  where 
the  electron  density  in  the  boundary  layer  plasma  will  be  lower  than 
anywhere  else  in  the  layer.  For  the  operating  frequency  of  Fig.  33, 

2  GHz,  the  electron  (tensity  which  would  result  in  a  plasma  frequency 
equal  to  the  operating  frequency  is  5  x  1010  electrons/cm3.  For  Fig. 

34  the  operating  frequency  is  40  GHz  and  the  electron  density  yielding 
a  plasma  frequency  equal  to  the  operating  frequency  is  2  x  1013  elec¬ 
trons/cm3.  In  both  Figs.  33  and  34  the  conversion  efficiency  decreases 
for  electron  densities  larger  than  that  electron  density  which  deter¬ 
mines  a  plasma  frequency  equal  to  the  operating  frequency. 

The  acoustic-to-electromagnetic  conversion  efficiency  is  plotted 
as  a  function  of  plasma  temperature  in  Fig.  35.  The  conversion  effi¬ 
ciency  (the  ratio  |^zt/p^|)  is  independent  of  the  temperature.  Even 
though  the  ratio  is  constant  for  varying  temperatures,  p^  will  be 
smaller  for  lower  temperatures;  hence,  the  effect  will  be  a  reduction 
in  the  transmitted  magnetic  field  at  lower  temperatures. 
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The  acoustic-to-electrouagnetic  conversion  efficiency  is  plotted 
as  a  function  of  the  dinensionless  parameter  in  Figr.  36,  37,  and 

38  for  three  different  electron  densities:  1011  cm-1,  iO,J  em'\  and 
10,:*  cn”3.  There  is  a  maximum  in  the  conversion  efficiency  at  a  fre¬ 
quency  slightly  above  the  plasma  frequency  and  the  efficiency  decreases 
continuously  at  higher  and  lower  frequencies.  It  is  more  difficult  to 
explain  the  trends  of  the  frequency  dependence  of  the  acoustic-to- 
electromagnetic  conversion  efficiency  by  physical  reasoning  than  to 
explain  the  frequency  dependence  of  the  electromagnetie-to-acoustic 
conversion  efficiency.  In  the  latter  case,  the  electromagnetic  wave 
actually  interacts  with  the  plasma  in  order  to  excite  organized, 
longitudinal  oscillations  of  the  eleccron  gas  as  a  medium .  However,  in 
the  case  of  an  acoustic  wave  exciting  a  transverse  electromagnetic 
wave  across  a  dielectric-plasma  boundary,  the  interaction  is  more 
difficult  to  visualize  since  the  organized  oscillations  do  not  exist 
in  the  dielectric.  Since  the  conversion  of  a  pressure  wave  into  an 
electromagnetic  wave  at  a  dielectric-plasma  interface  is  the  reverse 
process  of  the  conversion,  of  an  electromagnetic  wave  into  a  pressure 
wave  at  the  interface,  it  would  be  reasonable  to  expect  the  frequency 
dependencies  of  the  conversion  efficiencies  to  be  roughly  inverse. 

This  is  true  at  frequencies  below  the  plasma  frequencies  where  electro¬ 
magnet!  c-to-aco  us  tic  conversion  efficiencies  increase  and  acoustic-to- 
electromagnetic  conversion  efficiencies  decrease.  Also,  the  electro- 
magnetic-to-acoustic  conversion  efficiencies  have  a  minimum  for  w  o.  , 
whereas  the  acoustic-to-electromagnetic  conversion  efficiencies  have 
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Figure  37.  Acoustic-to~Electronagnetic  Conversion  Efficiency 
versus  (w/ui  )  „  (n  ■  10 1 3  cn“3) 
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a  Maximus  for  w  2  m  .  A  significant  result  is  the  fact  that  both  con- 
P 

version  efficiencies  decrease  for  higher  frequencies. 

The  acoustic-to-electromagnetic  conversion  efficiency  is  plotted 
as  a  function  of  the  dimensionless  parameter  v/tu  in  Fig.  39.  The 
effect  of  the  collision  frequency  on  the  efficiency  is  to  reduce  the 
magnitude  of  the  efficiency.  From  Eq.  (145),  if  v  +  »,  the  conversion 
efficiency  would  be  zero  (note  A  *  jw  +  v  in  Eq.  (145)).  me  effect 
of  the  collision  frequency  upon  the  acoustic-to-electromacnctic  con¬ 
version  efficiency  is  even  more  severe  than  upon  the  electroiia^netic- 
to-acoustic  conversion  efficiency  (Fig.  31).  The  effect  of  the  colli¬ 
sion  frequency  upon  the  acoustic-to-electromagnetic  conversion  effi¬ 
ciency  begins  for  collision  frequencies  well  below  the  operating 
frequency. 

Power  Conversion 

It  was  shown  in  Chapter  IV  that  the  amplitudes  of  any  traveling 
wave  pressure  solution  (or  superposition  of  traveling  waves)  of  the 

i 

coupled  wave  equations  were  negligible  compared  to  the  "forced  solu¬ 
tion,"  which  is  dependent  mainly  upon  the  coupling  term  on  the  right- 
hand  side  of  Eq.  (62).  It  is  these  traveling  (or  propagating)  pressure 
waves  with  wavelengths  of  the  order  of  magnitude  of  one  thousandth  of 
a  millimeter  which  are  required  to  obtain  the  resolution  needed  in  the 
very  short  (1  mm)  distance  from  the  wall  to  the  peak  electron  density. 
It  is  these  waves  for  which  the  conversion  efficiencies  have  been 
calculated  in  the  last  two  sections.  Even  though  it  has  been  shown 
that  the  amplitudes  of  these  waves  are  negligible  compared  to  the 


Figure  39.  Acoustic-tc-Electromagnetic  Conversion  Efficiency  versus  (v/w) 


"forced  solution*"  it  is  still  desirable  to  obtain  sow*  e*tirsate  of 
the  order  of  magnitude  of  the  amplitudes  of  these  waves.  It  will  ho 
shown  that  the  amplitudes  arc  functions  of  the  magnetic  field  strength 
however,  the  amplitudes  are  still  small  and  the  waves  axe  severely 
damped. 

Using  the  same  assumption  about  the  equation  coefficients  and 
g  (g  »  FH  ) ,  i.e.,  that  they  are  all  constants*  an  in  the  section 
Solutions  of  the  Coupled  Wave  Equations  in  Chapter  IV,  Che  complete 
solution  for  p  is  than  given  by 

V  >?*  u 

p  -  A'  e  +  B’  e  +  K/E  (146) 

where 

Yj  *  a  +  jB  M  i</E 

and 


Y2  =  -a  -  jB  -  , 

(ot  is  the  attenuation  constant  and  6  is  the  propagation  constant)  ar,d 
A'  and  B'  are  determined  from  boundary  conditions.  The  amplitude  of 
the  pressure  wave  traveling  in  the  positive  x  direction  is  B'  and  the 
amplitude  of  the  reflected  pressure  wave  (the  wave  traveling  in  the 
negative  x  direction)  is  A'. 

The  complex  amplitudes  A'  and  B’  will  be  calculated  using  the 
two  acoustic  boundary  conditions  from  Eqs.  (88)  and  (90)  for  a  plasma 
model  like  that  pictured  in  Fig.  5(b),  It  is  believed  that  this  model 
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Being  Eq.  (89)  tai  aoO£§  that  y  ?  v,  (since  x£  i#  nseb  eaaller 

P*  f 

tb»  tit*  other  tttm  to  iq.  (!K)))  gl 


V 


_  t  y~i  k.4 

-jA*  *J  *  y^*’  *  *  *  -T,*r  e  *  -  7,8’  *  x 


* 


where 


< 


n*s  been  e&eueed  to  be  aero  and  the  coapiete  pressure  solution  is  used 
iu  applying  the  botnriaxy  conditions.  Proa  the  preceding  equation 


(since  y2  •  -Y,) 


II? 


-wl 


(148) 


lit  leg  the  pi  as**  |h<raKt9iB  given  above,  th*  value  or 


<*/£> I 


ls»d 


froK  fig  (1?)  and  ncudcig  a  *  10  4  a"*1,  Eq«.  (14?)  and  (148)  btcaw 

-  i'  5  5  *  m~b 


m4 


*A)  ae  2.5  x  *Tn  . 

Stoat  I  A*  f  -  5  *  UT*\ 

18*1?  5  x  13”5  . 

Therefore ,  for  thia  profile,  the  amplitude  of  the  pressure  wave  propa¬ 
gating  mtmy  fro*  the  botndciy  is  5  x  10*s  newtons /a2  at  die  boundary, 
this  value  is  larger  than  die  "forced  solutions"  at  the  boundary; 
however,  die  propagation  constant  contains  at<  attenuation  factor  which 
causes  this  wave  to  be  dumped  out  very  rapidly  in  exactly  the  same 
manner  as  for  a  homogeneous  plasma  (Fig.  22).  The  value  of  *  is  a tore 
likely  to  be  105nfl  rather  than  lO4®”*,  in  which  case  { A* j  »  0  and 
the  wave  could  be  damped  out  in  about  one  hundredth  of  a  millimeter. 
For  a  =  104  m-1,  the  amplitude  of  the  pressure  wave  traveling  in  the 
negative  x  direction  is  2.5  x  10”11  newtons /m2  (at  x  »  0),  It  is  this 
amplitude  which  would  be  used  in  calculating  the  conversion  from  a 


fnuwe  vara*  into  m  elcctraMpetfc  v*w  at  to*  boundary. 


It  iffcan  tfaai  die 


wlatiea  (the  wywretial  or 


pnptt*tio|  vara*  factor)  for  the  ^kmkc  field  way  be  larger  (oear 


ft*  boaiaijr)  than  the  *(rt«4  wMUa*  be  co  toe  bo— <arf  c— ditto— . 
This  «w  t fee  cave  for  the  a— deal  solatl— s  also,  la  the  *c radical 


r,  the 


saiatloM,  the 


field  bad  a  "spika"  soar  *  -  0; 


dele&l—  ^jddly  converged  (ia  a  t**f  fitort  distance)  to  the  '‘forced 


eolation" 


the  affect  of  the  exponential  factor*  of  the  hr 


aoitticc  v—  sat  observable.  This  la  pres— ed  to  he  caved  b/ 


da— left  expe-lenced  bp  the  propagating  aodee.  The  Vpikes 


my  of  tor  f^enw  ill— trating  pleasure  solatia— 


(Figs.  13  to  It  — d  Figs,  20  — d  21) ,  chare  oelj  toe  "forced  eolutlooa* 
are  plotted. 

Ihe  power  la  —  elect roeaf— tic  vase  excited  in  toe  dielectric 
bp  a  pressure  v awe  incident  ft—  toe  pi— * a  will  a—  be  calculated. 

The  Pointing  vector  will  be  calculated  et  toe  bojadaip  and  toe  closed 
surface  shown  in  Fig.  W  will  be  —  ed  to  detendne  the  power  radiated 
in  toe  firlte  two-dt— nsional  geo-try.  The  width  (dx)  of  toe  volu— 
will  be  considered  to  approach  zero.  The  surface  area  of  toe  volu¬ 
te  toe  y-z  plane  la  S  —tar2. 

From  Pointing's  theorem  the  tl— -averaged  power  is 


*K}-* 

Closed  surface 


(149) 


where  ds  ■  -i  dydz  with  i  representing  the  unit  vector  in  the  x  direc¬ 
tion  and  the  asterisk  implies  a  complex  conjugate  operation.  Figure  32 


,  •.  whittsui  >,  *>  wh>  w'»vSr*Vu>.\4!u\vV\\' .  »,\\  *A  * . 
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iUwtfMct  tM  iftddtst  presses*  w*«  n£  ftt  excited  ehctrmpttlc 
v«w,  S{,  in  *i rr.  dieienric. 

For  •  eettici'.lj  poi«ri»d  tiectna^aetic  »m  ,  the  electric 
fielJ  in  the  diei  cirlc,  E^,  bm  ccmpmtmta  iM  Ey  Md  8f  bee  only 
a  z  csapeceac  B,.  For  Che  y  dependence,  of  If.  (140), 

2t  r  -  +  iEyS*  -  JExM*  ,  030) 

* 

where  E  .  E  «  and  H  ere  functions  of  x  only  since  the  e  7  ten  will 

jf  X  £ 

cancel  is  the  operation.  Substituting  Eq.  (ISO)  into  Eq.  (149)  gl 


F*  • 


since  E  and  B  axe  functions  of  x  only. 

y  * 


P  -  -  Re  |  -EyH*  |  S  (watts)  . 


(151) 


The  electric  field,  Ey,  in  the  dielectric  is  given  by 

.  3H 
1  z 


y  j<**d  3*  * 


or 


H  Y  .  Y,  X 

zt  dx  dx 

jU)C.  ® 


where  H  is  the  amplitude  of  the  excited  magnetic  field  strength  at 

Z  l 

x  ■  0.  Using  Eq.  (135)  fer  y^  and  evaluating  Ev  at  x  ■  0”  (just 


inside  the  dielectric)  gives 


E  -  -H  n,  cos  6, 
y  zt  c.  1 


(152) 
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Khar* 


A 

is  A*  characteristic  iapehoa 

StkaUtotlac  (152)  lute  Eq.  (251)  |i«M 

F  *  I  [8«  r,d  <*•  *lE*]  *  *  *W3> 

*  -  * 

ffSloitM  at  *  *  0  is  B^.  Squstiec  (153)  dsn  become 

'<*«*«]  *  • 

ft  2 

Sine*  HEtElt  i«  jvt  jSJ{l  which  is  a  ml  ^cwtity,  d«  ti*e- averaged 
power  is 

,2  2 

Fr vm  f««.  (145)  fits  ratio  tj  Ai?jJ  esn  be  doteraiced.  This  ratio 
will  be  represented  by  C2: 


•  *>7t  E. 

o  9  c 


of  tbs  dielectric. 


Then  P  can  be  expressed  as 

—  1  iii2 

P  »  y  nd  cos  6,C*  S  iPji  (watts) 


(154) 


Using  Eq.  (154) ,  the  power  In  an  electromagnetic  wave  excited  by 
an  incident  pressure  wave  will  be  calculated  for  the  following  param¬ 
eters: 

(1)  Angle  of  incidence  30° 

(2)  Operating  frequency 


30  GHz 
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(3J 

kd 

5.S 

(4) 

S  (1  i*2) 

4.3  x  10*4,  neter* 

(5) 

Plants  density  profile  ft) 

10 1  »  .  JO***  eg*3 

(b) 

Flueae  temperature  profile  ft) 

2000-5000*8 

'7} 

Collision  frequency 

4  x  10* 1  sac'1 

(S) 

Distance  to  the  peak  elect rca 
density 

1  *o 

As  stored  at  the  beglnalag  of  this  mti«k  tb*  appropriate  anplitode 
to  bo  ur«d  In  aMMtte-to^lKtiraugMtic  tosiettlai  colcnUtiou  U 
A*  r  tie  amplitude  of  the  pmacrs  vac*  trawling  in  the  at|ttlwt  a 
direction  «u»0.  For  cho  fnwttw  (Ivon  (boo*,  the  Mtlaat*  for 
A*  gives  by  Eq.  <148)  v«  fond  to  bo  2.$  x  30'* 1  newt nos /a2.  tta, 
|plj  *  *  l°~n  aowtoosAo2.  2bo  acoue tlc-toc-olacttwe gni tic  cm- 

verelen  efficiency,  C  In  Eq.  (154),  cat  be  estimated  fra*  Figs*  33 
and  34.  In  Pit.  33,  C  la  2.5  for  a  frequency  of  2  a*  aid  is  Pig.  34, 
C  is  1.5  for  a  frequency  of  40  GBt,  ubei*  an  electron  daoaity  of  Id** 
ci’J  was  used  to  determine  C.  For  an  operating  jt.vqcancy  of  30  GHx 
C  will  be  somewhere  batween  1.5  and  2.5,  so  lot  C  be  equal  Co  2  a£ 
an  approximation.  All  of  the  factors  In  Eq.  (154)  have  bean  deter¬ 
mined  and  the  ties- averaged  steady-state  power  can  be  calculated  as 
follcws : 

P  -  (0. 5) (202) (0. 866) (2)2 (6. 3  x  HT4)  <2.5  x  30'11)2 


which  gives 


P  -  1.4  x  10”22 


watts 


13»i*  i «  cn  ntnwly  mil  ptaw r  lcwi  ed  it  would  be  wry  bard  to 
distinguish  between  this  signal  sod  background  noise  due  to  spectral 
Mission  in  the  plans  layer,  etc.  Is  tame  of  power  relative  to 
1  milliwatt  (or  dBm)  the  power  level  is 

P  -  -188  dBm  , 

which  would  require  4  very  sensitive  receiver  indeed  for  detection. 

It  must  be  emphasized  that  the  attenuation  factor  was  assisted  to  be 
ID4  af1.  If  *he  attenuation  factor  is  actually  an  order  of  magnitude 
larger  than  tills  (which  it  more  than  likely  is),  the  power  given  by 
Eq*  (154)  would  be  virtually  zero. 


CK4P7ES  VI 


CUKCUSIOB 

The  propagiticrt  of  electron  acoustic  waves  in  an  inhomogeneous , 
warm,  lossy  plasma  layer  has  been  investigated.  The  coupled  wave 
equations  (Eqs.  (6l)  and  (62))  have  been  solved  by  a  finite 
difference  numerical  technique  in  order  to  determine  the  degree  of 
coupling  between  the  electromagnetic  wave  and  the  electron  acoustJ-. 
wave  in  the  inhof />geneous  region  of  a  warm ,  lossy  plasma  layer. 

For  the  plasma  profiles  and  operating  frequencies  considered,  it  was 
found  that  the  pre-'  jure  (or  acoustic)  solutions  of  the  coupled  wave 
equations  consisted  of  two  parts:  a  "forced  solution"  with  an  x- 
voriation  similar  to  the  coupling  coefficient  in  the  pressure 
equation  (Bq.  (62}),  and  a  "hoacgeneous  solution"  which  character¬ 
ized  the  propagating  electron  acoustic  wave.  It  is  shown  that  both 
parts  of  the  solution  are  directly  dependent  upon  (or  coupled  to) 
the  magnetic  field  in  the  inhomogeneous  plasma  (see  Iqs.  (123), 

(lho),  (1&7)  and  (1^6)}.  It  is  also  shown  that  for  the  plasma 
profiles  and  operating  frequencies  considered,  the  homogeneous  solu¬ 
tion  is  negligible  compered  to  the  forced  solution  except  possibly 
near  the  dielectric-p] tsma  interface.  Thus,  even  though  the  electron 
acoustic  wave  is  indeed  coupled  to  the  electromagnetic  wave,  the 
electron-acoustic  wave  i&  severly  damped  for  propagation  in  either 
direction  in  the  inhomogeneous  plasma.  In  the  approximate  analysis 
it  is  shown  that  the  "propagation  constant"  for  the  electron  acoustic 
wave  in  the  inhomogeneous  region  is  the  same  as  that  for  a  homogeneous 


12U 


m 

icsty  plassg*  (Sq,  {13C>)>  acd  the  ansj-iag  effect  is  due  to  tse  iar^f 
«  .Ilia  '.oo  frequency  cf  the  plassb 

It  is  also  shows  la  the  approximate  an&Iysi*  that  the  .teaiy- 
siate  t  ime«aversq,»<l  power  produced  by  the  cpefersioa  ©f  «  ore«*ure 
ware  let©  an  »Jectro«*gflj?tfe  wav*  at  &  dielectric  boundary  is  of  the 
order  of  1C-*2  watts  for  a  typical  plflsfia  profile  rod  operating 
frequency.  This  rasvlt  is  baaed  on  a  probe  area  cf  1.0  ir/ ,  as 
incident  electromagnetic  power  of  10  aw  end  an  operating  frequency 
of  30  GKx.  This  received  powsr  level  is  loo  lew  fer  reliable 
detection. 

In  addition,  it  is  shown  that  the  conversion  efficiently  for 
pressure  wave- c o-s Xec  t rooagn et  1  e  vavo  at  a  dielectric-plas** 
boundary  is  reduced  for  collision  frequencies  well  below  the  opera¬ 
ting  frequency  {see  Fig,  31).  High  collision  frequencies  will, 
therefore,  reduce  the  detection  capability  of  an  acoustic  probe 
system  by  damping  the  propagating  acoustic  wave  ar.d  by  reducing  the 
conversion  efficiency.  The  results  of  this  investigation  show  that 
even  for  an  inhomogeneous  plasma*  the  damping  effects  of  high 
collision  frequency  on  electron-acoustic  waves  is  drastic. 

For  the  plasma  profiles  considered,  and  using  a  steady-state 
plane  wave  analysis,  the  effect  of  the  acoustic  solution  upon  the 
electromagnetic  solution  is  shown  to  be  negligible.  This  is  true 
even  for  a  "warm”  plasma  layer  with  temperatures  as  hig*.  as  5«XX)CK. 
Therefore,  for  the  plasma  conditions  studied,  the  approximate 
solution  for  the  pressure  given  by  Eq.  (115),  used  together  with 
Eq.  (6i)  for  the  magnetic  field,  adequately  treats  the  effects  of 


coepressibi -ity.  ?his  is  valid  eisr-e  the  approximate  solution  for 
P  is  very  good  for  the  plasma  parameters  considered.  Thus  the 
effect  of  the  acoustic  portion  of  the  fields  in  the  plesm&  is  ta>ten 
care  of  by  substituting  p  from  the  approximate  solution  into  the 
right  hand  side  of  £q.  (6l).  In  this  case*  the  effect  on  the 
electromagnetic  fields  in  the  plasma  will  be  determined  by  solving 
Eq.  (cl 5  rather  than  the  tvc  coupled  second  order  differential 
equations . 
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APPENDIX  A:  APPROXIMATIONS  INVOLVED  IN  THE  APPROXIMATE  SOLUTION 

The  approximations  used  in  obtaining  the  approximate  solution  will 
be  investigated  here  for  the  following  parameters: 

n  :  1011  -  1014  cm*3 

o 

T:  3000  -  4000°K 
v:  6.5  x  10 J 1  sec*1 
f:  2  GHz 

0  j :  20° 

The  coefficients  for  the  differential  equations  were  computed  and 
their  maximum  or  minlmtxn  absolute  values  are  listed  below: 

|A|  6.0  x  102 

|fi|  <_  2.0  x  104 
|C|  <  1.5  x  105 
!d|  <_  6.0  x  105 
|E|  <_  1.8  X  1012 
1 F |  <  2.0  x  105 

The  characteristic  length  of  the  plasma  inhomogeneity  is  just  the  width 
which  was  2mm.  The  approximation  will  improve  as  d  is  decreased. 

Using  Eqs.  (108)  and  (109)  gives 

(2.5  x  105  +  3  x  105  +  2  x  104)  H  *  1.5  x  105  p 

z 

and 

(2.5  x  105  +  3  x  108  +  1.8  x  1012)  p  -  2  x  10fc  H 

z 

From  the  first  equation,  it  is  observed  that  the  first  two  terms  in 
the  parentheses  are  much  larger  than  the  last.  In  the  second  equation 
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the  third  tens  in  the  parentheses*  completely  dominates  the  first  two. 

In  the  approximate  solution,  only  the  dead n ant  terms  are  retained. 

It  is  now  necessary  to  compare  CF/E  with  1/d2  mad  A/d.  For  the 
profile  and  operating  frequency,  |e|  minimum  •  8*  10*°  and  the  maxi¬ 
mum  absolute  values  for  C  and  F  have  been  given. 

CF 
E 


i-  -  2.50  x  10 5 
d2 


-  3.75  x  102 

'max. 


-  4.00  x  10b 


Therefore  CF/E  Is  much  smaller  than  the  ether  two  terms  and  is  neglected 
in  the  approximate  solution.  In  the  approximate  fore,  the  equations 
then  becooe 


( 


«  0 


and 


Some  general  remarks  can  be  made  about  the  validity  of  the 
approximate  solution.  The  approximation  la  move  accurate  at  lower 
frequencies.  The  accuracy  la  also  improved  as  the  electron  density 
is  increased  or  die  plasma  thickness  is  decreased. 


APPENDIX  B:  COEFFICIENTS  OF  THE  COUPLED  HAVE  EQUATIONS 
FOB  A  PARTICULAR  CASE 
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The  real  and  imaginary  parts  of  the  coefficients  of  the  coupled 
vase  equations  ere  given  in  Table  II  for  the  following  plasma  profile: 


n  :  10s2  -  101J  cm"1 

o 

T:  2000  -  3500 *K 
d:  1.0  mm 
v:  101*  sec*1 
6,:  30* 

is  26  GH* 


the  coefficients  chao^e,  of  coarse,  for  different  plasma  profiles, 
angles  of  Incidence  and  operating  frequencies. 


NOTE:  R(A)  means  the  real  part  of  coefficient  A  and  1(A)  means  the 

imaginary  part;  the  same  notation  applies  for  the  other  coeffi¬ 
cients. 
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APPENDIX  C:  APPROXIMATION  FOR  THE  DERIVATIVE  OF  B 

z 

AT  THE  FAR  BOUNDARY 


The  expression  for  the  derivative  of  H  at  the  outer  boundary  for 

z 

the  plasma- free  space  condition  is  given  by 


dH 
_ z 

dx 


/eK  \ 

-  i 

x*d 

W  -  y«  ».d) 

where  p(d)  and  are  the  values  of  p  and  at  the  boundary.  It  is 
necessary  now  to  compare  the  two  terms  on  the  right-hand  side  of  the 
equation  above. 

£  P«)  <  yo>  »Id  • 


Substituting  typical  values  for  the  variables  above  (and  choosing  them 
so  that  the  left-hand  aide  would  be  as  large  as  possible  and  the  right- 
hand  side  would  be  as  small  as  possible)  gives 


10~19  103 
10 11  10-30 


?  3 

p(d>  <  (0.1) (10) 


10  pd)  <  . 

From  the  section  Solutions  of  the  Coupled  Wave  Equations  it  is  seen 
that  p(d)  is  always  at  least  four  orders  of  magnitude  below  Hz<j,  there¬ 
fore  the  inequality  above  holds  and  we  have 
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cK 


V«  '« “.d 


* 


and  to  a  very  good  approximation. 


xd  ' 


dU 
_ z 

dx 


-je  (d)  K  H 
J  p  ox  : 
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APfLbDlX  D:  ACOUSTIC  WAVELENGTHS  COME  AXED  TO  THE  DEBTE 
LENGTH  AID  THE  IXTEX-PAPTICLE  SPACING 


The  acoustic  warns  length  in  the  inhomogeneous  plana  will  be 
ccaputtd  at  the  equivalent  wavelength  is  a  lossy,  homogeneous  plena 
with  the  »mm  electron  density  aid  fpttatwe.  This,  of  course,  1* 
only  an  approximation  since  a  definite  wavelength  cannot  be  found  for 
a  pluses  wave  In  an  inhoaogeneous  medium  The  wave  length  la  computed 
as  follows: 


P* 


I 771 

VT  'i 


M2 

4-  _E  + 

u2 

o 


.  •*!  v 

u2  " 
o 


where  YpI  consist*  of  a  real  part  and  an  1  nag! nary  part  a>d  Is  just 
the  x  component  of  the  conplex  propagation  constant  of  the  plasaa  wave 
la  a  homogeneous  plasaa.  The  plaaea  wavelength  In  the  x  direction  la 

v2*4{v}  • 


where  In  Indicates  that  the  Imaginary  part  of  y  la  to  be  used. 

P* 

The  puzpoee  of  computing  a  plasma  wavelength  is  to  determine 
if  the  plasma  wavelength  la  larger  or  smaller  than  the  Debye  length 
in  the  plasma.  The  Debye  length  la  computed  from  the  following 
expression: 


The  placw  vawJengtb  is  also  compared  to  the  inter-particle 


1JI 


spacing  which  la  given  by  n^" . 

The  three  q tun tl ties  described  above  (plmna  wavelength,  Debye 
length,  and  inter-particle  spacing)  were  computed  an  fanctinna  of  the 
distance  Into  the  piaana  layer  in  ?igs.  *1  to  44  for  the  four  plana 
profiles  of  interest. 

it  is  noted  that  in  Figs.  42  end  4i,  the  plana  wareletgtfa  la 
saaller  than  die  Debye  length  in  a  very  nartoe  region  close  to  the 
wall.  This  is  due  to  the  low  electron  densities  In  this  region.  1c 
Fig.  43  it  is  observed  that  a  decrease  in  operating  frequency  can 
increase  the  plana  wavelength  to  s  raise  larger  then  the  Debye 
length.  The  plana  wavelength  is  larger  than  tne  inter-particle 
spacing  in  all  cases. 


LCMTH  SOLI  tot tt) 


Figure  41.  Acoustic  Wavelength,  Debye  Length  and  luter-Particle 
Spacing  (n  :  1012  to  1013  c«r3) 


Figure  44.  Acoustic  Wavelength,  Debye  Length  and  Inter-Particle 
Spacing  (n  :  1013  to  1014  cm-3) 
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AJPE51IX  E:  riGXT.AL  '"'WFUTEF  FFOGBA.V 

Calculations  of  the  electromagnetic  sad  pressure  fields  in 
the  inhomogeneous  pi  as  as  sheath  were  performed  osi  the  CLC  6600 
computer  at  the  Aiv  Force  Veepcns  Laboratory,  Ki it land  AFB,  lev 
Mexico.  Based  os  these  calculations  for  various  plasca  sheaths 
vitli  large  peak  electron  densities  and  high  collision  frequencies, 
it  was  eonciui ei  in  this  study  that  the  snail  wavelength  propa¬ 
gating  acoustic  waves  were  daaped  so  severely  that  reliable 
detection  of  the  reflected  acoustic  waves  would  be  impossible. 

For  less  dense  plascas  or  for  plasnaa  with  such  lover  collision 
frequencies,  reflected  acoustic  waves  nay  provide  valuable  diag¬ 
nostic  information.  The  computer  progrss  developed  for  this  work 
would  apply  equally  well  for  these  types  of  plasnas  and  the  program 
is  included  for  this  reason.  A  FOHTRAS  listing  of  the  progrss  and 
subroutines,  along  with  sample  input  and  output  data,  are  given  in 
this  appendix. 

Program  FSCDE1  uses  the  finite  difference  numerical  solution 
process  described  in  Chapter  IV  to  solve  the  two  equations: 


d2H  .  dH 

+  A  ~  +  BE 
.  o  dx 
ax* 


— £  +  D_p  +  E  =  FH 

dx2  d£  p 


where  A,  B,  C,  D,  E  and  F  are  functions  of  the  plasma  medium  and 
vary  with  distance  x.  The  program  is  divided  into  subroutines  which: 


(1)  calculate  the  boundary  conditions; 

(2)  calculate  the  coefficients; 

(3)  calculate  the  initial  "guesses"  for  H  and  p; 

(i*)  iterate  the  numerical  solution; 

(5)  set  the  output  data  format; 

(6)  calculate  the  reflection  coefficient. 

There  is  one  basic  restriction  on  the  program  as  it  is  presently 
written.  The  electron  density  profile  must  be  linearly  increasing 
from  a  value  at  the  wall  to  a  peak  value  at  3ome  distance  from  the 
wall.  Both  cf  the  models  in  Fig.  5  can  be  used  if  appropriate 
boundary  conditions  at  the  location  of  the  peak  electron  density  are 
used  in  the  subroutine  calculating  the  boundary  conditions. 

Listed  below  are  the  FORTRAH  nas.es  assigned  to  the  input  and 
output  variables  of  Program  USCBE1,  along  with  their  definitions  and 


units . 

INPUT 

FORTRAH  NAME 

DEFINITION 

UNITS 

AKD 

Relative  dielectric  constant  of 

Dimensionless 

the  dielectric  material. 

AH1 

Electron  density  at  the  wall. 

electrons /m 3 

AN2 

Peak  electron  density. 

electrons/m3 

ANU 

Collision  frequency. 

sec  1 

BLD 

Distance  from  the  wall  to  the 

meters 

peak  electron  density. 

DEGREE 

Angle  of  incidence  (0^). 

degrees 

FREQ 

Frequency  of  electromagnetic  wave. 

Hz 

FORTRAN  SAMS  DEFINITION  UR  ITS 

MAXIT  Xaxitsum  number  of  iterations-  integer 

5FT5  Number  of  points  on  the  x-axis  integer 

for  the  numerical  solution. 

TE1  Elect ron  temperature  at  the  wall.  °K 

TE2  Electron  temperature  at  x  *  BLD  °K 


Mete:  AN1,  AN2,  AMU,  TE1  and  TE2  are  defined  in  the  program 
under  PROFILE  DEFINITION.  The  variables  AKD,  BLD,  FREQ  and  DEGREE 
are  read  into  the  program  under  FORMAT  HE20.2.  The  variable  HPTS 
is  read  into  the  program  under  the  FORMAT  110. 

OUTPUT 

FORTRAN  NAME  DEFINITION  UNITS 

ABSH(I)  |H|,  magnitude  of  magnetic  intensity  amperes/meter 

intensity  at  point  I. 

ABSF(l)  jp|,  magnitude  of  pressure  field  at  newtons /meter2 

point  I. 

H(l)  Magnetic  intensity  at  point  I,  a  anrpere/meter 

complex  variable. 

I  Point  along  x-axis,  1=1,  2,...NPTS.  integer 

P(I)  Pressure  field  amplitude  at  point  I,  newtons /meter2 

a  complex  variable. 

XA(l)  Distance  value,  x,  at  point  I.  meters 

The  program  will  also  print  out  the  coefficients  o^  the 
coupled  differential  equations  as  functions  of  I.  In  this  printout, 
R(A)  represents  the  real  part  of  coefficient  A  and  l(A)  its  imaginary 
part.  A  similar  notation  is  used  for  the  other  coefficients.  In  the 


output  listing,  x  is  XA(I),  R(H)  is  the  real  part  cf  H(l),  1(H) 
is  the  imaginary  part  of  H(I),  H(P)  is  the  real  part  of  P(l),  I ( P ) 
is  the  imaginary  part  of  P(I),  ABSH  is  the  absolute  value  of  H(I) 
and  ABSP  is  the  absolute  value  of  P(I). 
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